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20.  ABSTRACT  (Continued) 

nitrate  (PETN)  explosives  were  used.  These  various  experiments  incorporated 
large  differences  in  explosive  specific  energy  soil  cohesiopjr,'  soil 
angle  of  internal  , friction  ♦  ,>  charge  weight  -Vr and~gravitational  acceleration, 
Experimental  results  show  substantial  differences  in  cratering 


phenomena. 
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The  gravity-scaled  yield  [  ^  =  (g/Qe) (b/6)  ].  previously  used  to  correlate 

centrifuge  cratering  results  in  dry  Ottawa  sand,  has  been  modified  to 
accurately  account  for  the  effects  of  variations  in  c  and  ♦.  When  cratering 
efficiency  is  plotted  as  a  function  of  this  strength-gravity-scaled  yield 
TU,all  experimental  results  fall  on  a  single  curve. 


These  results  show  that~Jt  small  actual  yields  soil  cohesion  dominates  the 
cratering  process.  At  iTrger  yields,  the  cohesion  has  little  effect,  but 
the  angle  of  internal  friction  becomes  an  important  variable  because  of  the 
large  overburden.  The  behavior  of  a  material  such  as  alluvium,  having 
nonzero  values  for  both  and  <$  ,  shows  cube-root  scaling  at  low  yields  but 
becomes  asymptotic  to  dry  sand  behavior  at  large  yields.  Neither  behaves  as 
a  zero-strength  material,  and  neither  shows  quarter-root  scaling  for  scaled 
yields  approaching  megatons  of  TNT.  ,  / 
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SUMMARY 


The  use  of  increased  gravity  as  a  necessary  and  sufficient  condition 
to  correctly  simulate  large-scale  cratering  events  with  small  laboratory 
experiments  was  theoretically  and  experimentally  verified  in  the  previous 
year's  work.  During  that  program  cratering  experiments  were  carried  out  for  a 
range  of  almost  11  decades  of  gravity-scaled  energy  in  dry  Ottawa  sand.  In 
addition,  a  highly  successful  simulation  of  the  JOHN  I E  BOY  nuclear  event  was 
performed  by  a  combination  of  a  nuclear  high-explosive  ( Ht£ )  equivalence 
generated  by  a  code  calculation,  together  with  the  small-scale  simulation  of 
the  resulting  HE  event  on  the  centrifuge. 

The  program  reported  here  was  intended  to  investigate  variations  in 
cratering  phenomena  due  to  differing  cratering  media.  In  particular,  parts  of 
the  program  included: 

1)  a  series  of  experiments  in  a  nominal  4-percent  moisture  desert 
alluvium  over  the  maximum  gravity-scaled  energy  range  of  about 

11  decades, 

2)  experiments  in  alluvium  with  moisture  content  from  near  zero  to 
7-1/2  percent, 

3)  experiments  in  a  fully  saturated,  dense  Ottawa  sand,  and 

4)  experiments  in  an  oil-base  modeling  clay. 

All  of  the  experiments  were  conducted  with  half-buried  spherical  charges. 

In  addition  to  the  experimental  proqram,  supporting  theoretical 
analyses  we^e  conducted  to  organize  and  interpret  the  experimental  results  and 
to  generate  a  rational  prediction  method  for  other  materials.  A  theory  was 
constructed  which  correlates  the  experimental  results  using  a  single  gravity- 
size-strength  parameter,  This  parameter  correctly  measures  the  equivalence 

of  all  combinations  of  gravity,  charge  size,  medium  and  charge  properties 
considered. 

The  theory  identifies  two  distinct  regimes.  For  a  given  soil  and 
charge  type,  a  cohesion-dominated  regime  exists  for  small  charge  size.  In  this 
regime,  cube-root  scaling  of  crater  volume  is  predicted.  The  value  of  the 
cratering  efficiency  in  this  regime  depends  directly  on  the  cohesion  of  the 
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material.  If  cohesion  is  zero,  such  as  for  dry  granular  materials,  this  regime 
is  not  observed. 

A  transition  to  a  1 ithostatic-pressure-dominated  regime  occurs  at  a 
value  of  charge  size  which  depends  upon  the  relation  between  the  cohesion  and 
the  angle  of  internal  friction.  For  charge  sizes  sufficiently  larger  than  this 
value,  the  cratering  efficiency  decreases  with  increasing  charge  size.  This 
decrease  approaches  a  straight  line  on  a  log-log  plot.  However,  a  decrease  to 
the  classical  “quarter-root"  scaling,  suggested  in  the  literature,  is  not 
observed.  This  failure  is  attributed  to  the  dependence  of  the  coupling 

i 

efficiency  on  the  specific  energy  of  the  explosive  and  on  the  other  variables 
included  in  the  theory. 

The  results  of  this  program  illustrate  the  utility  of  the  centrifuge  as 
a  tool  to  explore  the  dependence  of  cratering  phenomena  on  a  multitude  of 
parameters  in  a  precise  and  inexpensive  way.  The  results  of  the  desert  alluvium 
experiments  explain  the  discrepancy  between  the  small-scale  experiments  at  1  G 
and  the  large  field  shots,  and  furthermore  clarify  the  variations  found  in  1-G 
experiments  in  different  materials  such  as  saturated  sand. 
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SECTION  1 
INTRODUCTION 


Crater  size  and  shape  depends  upon  both  the  soil  type  and  the  explosive 
type  for  a  fixed  source  geometry.  Differences  are  attributed  to  the  influence 
of  material  properties.  Significant  variations  have  been  observed  in  the  many 
field  experiments  in  materials  ranging  from  water,  on  one  extreme,  to  granite, 
on  the  other,  and  for  energy  sources  consisting  of  various  nigh  explosives, 
nuclear  explosives,  and  hypervelocity-impact  projectiles. 

The  data  scatter  inherent  in  large-scale  field  events  has  made  it 
difficult  to  determine  quantitatively  how  medium  and  source  properties  affect 
cratering  mechanisms.  On  the  other  hand,  the  more  precise  laboratory  tests  at 
small  scale  do  not  always  reproduce  the  phenomenology  of  large  scale  events. 

Previous  results  of  a  theoretical  and  experimental  program  (Schmidt  and 
Holsapp1e,  1978a)  demonstrate  that  large  explosive  events  can  be  directly 
simulated  at  small  scale,  if  the  test  is  performed  at  elevated  gravity.  These 
conditions  are  obtained  by  performing  the  experiments  on  a  centrifuge,  with  a 
gravity  variation  from  10  to  greater  than  800  G.  Both  theoretical  and 
experimental  results  verify  that,  for  large-  and  small-scale  tests  performed  in 
the  same  material  with  the  same  explosive,  all  aspects  of  the  cratering 
phenomena  are  correctly  simulated.  A  possible  anomalous  effect  of  centrifuge 
testing,  e.g.,  Coriolis  force,  has  not  proven  to  be  significant.  The  energy  of 

3 

the  simulated  large-scale  test  is  equal  to  g  E,  where  E  is  the  energy  of  the 

small-scale  test  conducted  at  elevated  gravity  g.  Consequently,  a  variation  of 

gravity  from  1  to  500  G  allows  a  simulation  of  explosive  events  over  a  range  of 
3  6 

(500)  =  125  x  10  ,  or  more  than  8  decades  of  energy. 

Extending  this  earlier  work,  experiments  were  successfully  performed  in 
both  dry  and  wet  Ottawa  sand,  desert  alluvium  of  various  moisture  content,  and 
an  oil-base  clay  during  the  current  program.  Both  lead-azide  (PbN^)  and  PETN 
explosive  charges  were  used.  These  experiments  gave  data  over  a  large  range  of 
soil  strengths,  explosive  properties,  charge  mass  and  gravitational 
acceleration.  The  cratering  results  are  reproducible  and  consistent  with 
existing  field  events.  These  data,  together  with  a  consideration  of  the 
variations  in  material  properties  and  energy  sources,  provide  ar  empirical  base 
for  the  construction  of  a  theory  that  predicts  the  effects  of  media  strength  and 
energy-source  properties  on  crater  volume,  ra  ius  and  depth. 
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SECTION  2 

EXPERIMENTAL  RESULTS 


2-1  CENTRIFUGE  DESCRIPTION 

The  Boeing  600-G  geotechnical  centrifuge  was  used  in  this  study.  This 
machine  has  a  dynamic  load  rating  of  60.000  G-kg  (66  G-tons)  at  620  rpm  and  was 
constructed  using  the  aerodynamic  housing  and  main  shaft  assembly  from  a  Gyrex 
Model  2133  centrifuge.  The  rotor  was  designed  and  fabricated  by  the  Boeing 
Company  to  incorporate  symmetric  swing  baskets  for  geotechnic  applications.  The 
arm  radius  to  the  fully  extended  base  plate  is  139.7  cm.  The  maximum  payload 
mass  is  250  kg  on  each  rotor  end.  An  overall  view  is  shown  in  Fig.  2.1;  the 
details  of  the  swing  basket  and  soil-sample  container  are  shown  in  Fig.  2.2. 
Power  is  provided  by  a  30-horsepower  Eaton  Dynamatic  Model  ACM-326-910P  drive 
unit  incorporating  an  adjustable  speed,  constant -torque  eddy-current  clutch. 
The  unit  has  electrical  dynamic  braking,  allowing  shut-down  from  maximum  rpm  in 
less  than  30  seconds.  The  constant  speed  motor  and  variable  drive  unit  are  shock 
mounted  and  coupled  to  the  main  shaft  with  a  belt  to  minimize  vibration. 

The  rotor  shaft  is  equipped  with  24  slip  rings  For  instrumentation 
channels,  three  220-V. a .c .-power  slip  rings  and  a  hydraulic  slip  ring  which  can 
accommodate  either  gas  or  liquid.  A  pair  of  motor-driven  Nikon  F2  35-mm  still 
cameras  are  hub-mounted  in  a  stereo  configuration.  These  cameras  provide 
stereophoto  coverage  of  the  number-one  rotor  end  with  a  maximum  framing  rate  of 
six  per  second.  A  16 -mm  movie  camera  is  also  hub-mounted  and  operates  at  up  to 
500  frames  per  second,  giving  movie  coverage  of  the  number-one  rotor  end. 
Illumination  for  all  cameras  is  provided  by  three  600-watt  quartz-halogen  lamps. 
An  alternate  scheme  uses  Sylvania  FF-33  long-duration  ( 3- s*  c )  flash  lamps. 

2-2  RESULTS  OF  SMALL-SCALE  EXPERIMENTS 

Various  small-scale  explosive  experiments  were  performed  on  the 
centrifuge.  Spherical  explosive  charges  manufactured  by  the  R.  Stresau 
Laboratories  of  Spooner,  Wisconsin  were  used  exclusively.  These  included  three 
different  sizes  of  PETN  (0.49  qm,  1.34  gm  and  4.08  gm)  and  one  of  PbNg  (1.70 
gm).  All  were  centrally  initiated  by  applying  40  V.d.c.  to  a  notched 
0.0127-cm-diameter  tungsten  wi^e.  The  PETN  charges  contained  a  nominal  0.130-gm 
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;le  container  and  stereo-camera  configuration. 
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concentric  sphere  of  silver  azide  to  facilitate  ignition.  In  each  test  the 
explosive  sphere  was  half-buried  by  excavating  a  small  hole  after  screening  the 
final  surface  on  the  soil  sample. 

The  soils  tested  included  dry  Ottawa  sand,  desert  alluvium  obtained 
from  Kirtland  AFB  (with  moisture  content  from  near  zero  up  to  7.5  percent)  and 
an  oil-base  modeling  clay.  In  addition,  several  shots  were  made  in  a  nearly 
saturated  Ottawa  sand  demonstrating  the  feasibility  of  conducting  such 
experiments.  Details  of  the  experimental  technique  describing  soil  sample 
preparation,  determination  of  explosive  properties,  test  procedure  and  crater 
measurement  scheme  is  given  by  Schmidt  and  Holsapple  (1978a). 

A  complete  list  of  experimental  conditions  and  the  resulting  crater 
dimensions  is  given  in  Tables  2.1  through  2.4.  Appropriate  nondimens ional 
quantities  to  be  used  in  the  analysis  to  follow  are  also  tabulated.  The 
analysis  of  these  data  is  considered  in  the  following  sections. 

2-3  SATURATED-SAND  EXPERIMENTS 

The  final  series  of  shots  was  designed  to  assess  the  feasibility  of 
using  a  centrifuge  to  examine  cratering  behavior  in  saturated  media.  Four  shots 
in  Ottawa  sand  at  100  G  were  conducted  first.  The  first  two  (Run  26,  Table  2.4) 
provided  a  comparison  of  grain  size  effects  for  comparable  dense  sands.  Ottawa 
Sawing  sand  (nominal  range  of  grain  size  0.30-0.60  mm)  was  wetted  to  a  final 
wet  density  of  2.08  gm/cc.  This  was  contrasted  to  Ottawa  Banding  sand  (nominal 
range  of  grain  size  0.10-0.20  mm)  that  was  wetted  to  a  final  wet  density  of 
2.06  gm/cc.  Identical  charges  consisting  of  1.70  gm  spheres  of  lead  azide  were 
half  buried  in  each  sample.  The  crater  formed  in  the  coarser  grain  Sawing  sand 
was  2.23  times  larger  in  volume  than  that  produced  in  the  Banding  sand. 

A  Photo-Sonic  16 -mm  movie  camera  run  at  400  frames/sec  was  used  to 
record  the  shot  in  the  coarser  Sawing  sand.  The  camera  was  started  after  the 
centrifuge  was  a1;  the  appropriate  rpm  just  prior  to  firing,  providing  a  preshot 
reference.  The  centrifuge  was  held  at  100  G  for  60  seconds  after  firing  to 
observe  any  post-excavation  shape  modification,  due  to  possible  liquefaction  or 
slumping.  The  400-foot  film,  which  recorded  for  42  seconds  showed  no  evidence  of 
any  shape  change  following  crater  formation.  In  addition  to  the  movie  coverage, 
a  stereo-still-picture  sequence  confirms  that  no  significant  shape  change  took 
place  during  deceleration  of  the  centrifuge.  This  is  especially  significant 
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because  the  water  table  goes  through  a  shape  change  as  the  sample  is  accelerated 
from  1-G  conditions.  In  the  centrifuge  environment  the  water  conforms  to  a 
cylindrical  ceopotent i al  surface,  whereas  the  sand  surface  remains  flat. 
Therefore  to  achieve  a  water  table  the'  is  tangent  to  the  ground  surface  at  shot 
point  (127.1  cm  from  the  centrifuge  axis),  the  sand  surface  must  be  2.1  cm  below 
the  rim  of  the  46-cm-di ameter  container.  Under  these  conditions  the  water 
table,  initially  sufficiently  above  the  sand  surface,  will  contour  and  become 
tangent  to  the  surface  along  the  diameter  of  the  sample  that  is  parallel  to  the 
centrifuge  axis  of  rotation.  Excess  water  will  run  over  the  edge.  An 
alternative  to  this  configuration  would  be  to  contour  the  sand  surface  to 
conform  to  the  cylindrical  geopotential  surface.  In  either  case,  however,  there 
will  be  a  relative  motion  of  the  water  as  the  centrifuge  rotor  accelerates  from 
rest.  For  large  enough  arm  radius,  this  curvature  effect  can  be  neglected  as 
long  as  the  relative  water  motion  does  not  disturb  the  sand  upon  start  up  nr 
shut  down  of  the  centrifuge. 

The  second  pair  of  shots  (Run  27,  Table  2.4)  was  a  carefully  controlled 
reproducibility  test.  Ottawa  Flintshot  sand  was  pluviated  to  a  dry  density 
state  of  1.802  gm/cc.  Water  was  carefully  added  to  provide  a  final  wet  density 
of  2.113  gm/cc.  Within  the  accuracy  of  these  measurements  and  an  assumed  grain 
specific  gravity  of  2.65  gm/cc,  the  degree  of  saturation  was  greater  than  97 
percent.  The  goal  was  to  achieve  sand  states  comparable  to  those  of  Piekutowski 
(unpublished  data)  allowing  direct  comparison  of  the  centrifuge  data. 

The  results  indicate  a  high  degree  of  reproducibility.  The  crater 
radius  and  the  crater  depth  were  within  1  percent,  respectively,  for  the  two 
craters.  The  crater  volumes  differed  from  each  other  by  15  percent.  The 
coefficient  of  variation  (c.o.v.)  for  the  two  10U-G  centrifuge  shots  was  +11 
percent.  This  compares  favorably  with  c.o.v.  of  +1  3  percent  for  the  six  1-G 
shots  performed  by  Piekutowski.  Considering  the  complexity  of  charge  placement 
and  other  experimental  constraints,  associated  with  the  centrifuge,  these 
results  indicate  that  cratering  experiments  in  saturated  (or  nearly  saturated) 
media  can  be  satisfactorily  performed  at  high  G  in  the  centrifuge.  On  Run  27 
the  movie  camera  framing  rate  was  50  per  second  permitting  sufficient  recording 
time  to  follow  the  deceleration  phase,  confirming  the  observations  of  the 
previous  run  regarding  crater  stability. 

The  final  pair  of  shots  extended  the  range  of  gravity  variation  to  500 
G  and  provided  a  data  point  at  10  G.  Some  water  wave  washing  of  the  crater  was 
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observed  on  the  movie  coverage  of  the  1 0 -G  shot  (28-X).  This  filled  in  the 
crater  somewhat  creating  the  asymnetry  seen  in  the  profile  plots  in  Appendix  A. 

2-4  STRESS-WAVE  MEASUREMENTS 

A  series  of  soil-stress  measurements  were  made  using  carbon  gages 
placed  at  known  distances  from  the  explosive  charge.  Thirteen  shots  were  made 
under  various  operating  conditions.  Some  were  fired  with  the  centrifuge  not 
running.  These  provided  a  check  on  the  slip  rings  in  the  static  condition 
versus  direct  connection  to  the  oscilloscope,  bypassing  the  rings.  Six  stress 
gage  shots  were  made  with  the  centrifuge  in  operation  at  200  G  and  at  500  G. 
These  were  in  addition  to  the  crater  formation  shots  discussed  below. 

About  half  of  the  shots  were  successful  in  producing  stress  profiles. 
The  remaining  shots  failed  due  to  premature  trigger  problems  due  to  electrical 
coupling  in  the  slip  ring  circuitry  from  the  charge  firing  current.  A  3-mm 
diameter  piezo  pin  in  contact  with  the  charge  works  very  we i I  as  a  scope  trigger 
and  provides  a  reliable  time  reference  for  shock  initiation  into  the  soil.  From 
this  preliminary  series  of  tests,  it  is  concluded  that  dynamic  on-board 
measurements  can  be  made.  Appropriate  calibration  methods  need  to  be  developed. 
This  includes  not  only  stress  amplitude  but  a  technique  to  place  the  gage 
without  disturbing  the  sample  as  well  as  knowing  its  exact  location.  The  slip 
ring  transmission  appears  adequate  but  cross  talk  and  electrical  noise 
can  be  a  problem  under  some  circumstances.  This  can  be  eliminated  using 
on-board  recording  systems  which  are  currently  under  development. 

2-5  SOIL  MATERIAL  PROPERTY  TESTS 

To  facilitate  the  correlation  of  cratering  with  soil  material 
properties,  certain  laboratory  tests  were  performed.  The  bulk  of  the  tests  were 
performed  under  subcontract  to  Shannon  and  Wilson,  Inc.,  Geotechnic  Consultants, 
Seattle,  WA.  These  consisted  of  unconsol i dated-undrai ned  triaxial  tests  and 
also  some  direct  shear  tests.  The  reported  laboratory  test  data  supplied  by 
Shannon  and  Wilson  are  included  in  Appendix  R. 

Mechanical  property  data  were  obtained  for  Kirtland  Air  Force  Base 
(KAFB)  alluvium  at  various  moisture  contents,  dry  Ottawa  sand  of  three  different 
grain  sizes  and  Permoplast  oil-base  clay.  Two  series  of  triaxial  tests  were 
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performed,  the  first  was  conducted  at  high  confining  pressures,  up  to  27.5  bars. 
A  second  series  was  performed  at  lower  confining  pressures,  on  the  order  of  1.4 
bars  and  lower.  From  these  latter  tests,  soil  cohesion  and  angle  of  internal 
friction  were  determined  for  use  in  the  analytical  model  for  crater  volume 
(developed  in  Section  3),  incorporating  soil  material  properties.  The  direct 
shear  tests  were  performed  to  bracket  the  zero-confi ning- pressure  cohesion 
intercept.  These  values  are  considerably  lower  than  the  extrapolated  intercept 
from  the  low-pressure  triaxial  tests.  This  may  be  due  in  part  to  the  nonuniform 
stress  state  encountered  in  this  type  of  test.  For  calculational  purposes,  the 
direct  shear  values  were  assumed  to  be  a  lower  limit  for  the  cohesion  and  the 
low-pressure  triaxial  test  to  be  an  upper  bound.  The  hi gh-pressure  triaxial 
tests  show  some  reduction  of  tan  4,  which  is  commonly  observed  for  soils. 
Representative  values  for  all  three  soil  types  are  given  in  Sections  3  and  5. 
The  actual  test  results  are  included  in  Appendix  B. 
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SECTION  3 

A  MATERIAL  STRENGTH  MODEL  FOR  CRATER  VOLUME 


3-1  DIMENSIONLESS  PARAMETERS 

The  basis  for  direct  simulation  of  large  field  events  using  small 
laboratory  experiments  at  large  gravity  results  from  an  analysis  of  the 
similarity  requirements  between  the  two  experiments.  A  complete  similarity 
analysis  of  the  governing  equations  of  continuum  mechanics  in  their  general  form 
has  been  presented  by  Schmidt  and  Holsapple  (1978a,  1980)  and  will  not  be 
reproduced  here.  When  restricted  to  experiments  with  a  finite  set  of  governing 
variables,  a  set  of  dimensionless  parameters  can  be  derived.  Two  experiments 
are  similar  when  each  of  these  dimensionless  parameters  has  the  same  value  for 
both  experiments. 

For  the  application  here,  consider  the  case  of  a  half-buried 
(zero  depth  of  burst)  spherical  explosive  detonated  in  a  homogeneous  soil 
medium.  The  dependent  variable  of  interest  is  the  apparer  crater  volume  V.  It 
is  assumed  to  depend  on  soil  density  p,  soil  strength  Y,  explosive  mass  W, 
explosive  specific  energy  Qg,  explosive  mass  density  4 ,  and  gravity  g. 
Furthermore,  a  Mohr-Coulomb  strength  model  for  soil  is  assumed 

Y  =  c  +  P  tan  ♦  (3.1 ) 

so  that  the  strength  Y  under  confining  pressure  P  is  determined  by  the  cohesion 
c  and  the  angle  of  internal  friction  Therefore,  the  dependent  variable  V  is 
assumed  to  be  determined  by  seven  independent  variables,  p,  c,  ♦ ,  W,  Qg,  4,  g, 
for  a  total  of  eight  variables.  These  eight  variables  can  be  formed  into  five 
dimensionless  groups  defined  and  referred  to  as  follows: 


< 

(I 

IS 

(cratering  efficiency) 

(3.2a) 

a  W  1/3 

2  =  %  W 

(gravity-scaled  yield) 

(3.2b) 

p 

3  =  T 

(mass  density  ratio) 

(3.2c) 
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(cohesion  parameter) 


(3. 2d) 


*4'  bfa 

=  tan  4>  (internal  friction  parameter)  (3.2e) 

rherefore,  the  dependent  pi -group  must  be  determined  by  the  remaining  four, 

*y  =  F  (^2 *  n3»  1,4*  11 0 )  •  (3.3) 

For  a  given  soil-explosive  combination,  "j,  and  are  fixed,  and  *v  becomes 
a  function  solely  of  "o*  Two  experiments  in  the  same  material  using  the  same 
explosive  are  similar  whenever  ^as  the  same  value  and  they  will  have  the  same 
value  for  the  cratering  efficiency  *  .  These  two  experiments  can  differ  in 

v  3 

explosive  mass  W  and  in  gravity  g,  but,  in  order  to  be  similar,  the  product  g  W 
must  be  the  same.  Thus,  an  experiment  with  four  grams  of  explosive  at  500  G  is 
similar  to  and  directly  simulates  an  experiment  of  500  tons  of  the  same 
explosive  at  1  G. 

The  function  F  in  Eq.  3.3  can  be  determined  using  experimental  results 
and  its  dependence  upon  the  explosive  and  the  soil  properties  is  the  focus  of 
this  report.  The  differences  in  the  observed  ^  dependence  for  the  various 
soil-explosive  pairs  determine  the  influence  of  the  material  property  parameters 
and  The  form  of  this  latter  dependence  is  suggested  by  a 

consideration  of  the  fundamental  physics  of  the  phenomena,  as  will  be  shown. 

3-2  EFFECTS  OF  MATERIAL  STRENGTH  AND  EXPLOSIVE  PROPERTIES 

A  theory  that  unifies  the  experimental  results  presented  in  Fig.  3.1 
is  desired.  This  theory  must  account  for  and  quantify  the  variations  due  to 
medium  and  explosive  properties. 

The  previously  identified  variables  segregating  the  different 
materials  include  the  density  p,  cohesion  c  and  angle  of  internal  friction  $  of 
the  soil  medium,  and  the  density  6  and  specific  energy  Qg  of  the  explosive. 
These  variables  enter  into  the  three  dimensionless  pi-groups  and 

The  process  of  cratering  by  a  high  explosive  is  a  complex  result  of 
the  kinematics  and  dynamics  initiated  by  the  detonation.  This  detonation 
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process  is  commonly  modeled  as  a  shock  propagating  through  the  explosive,  with 
the  detonation  products  assumed  to  behave  as  a  perfect  gas.  The  pressure 
behind  the  shock  is  given  by  the  Chapman-Jouguet  pressure 

PCJ  •  2(y  -  l!«Qe  (3.4) 

where  y,  the  ratio  of  the  specific  heats,  is  taken  to  be  constant. 

The  effects  of  this  shock  on  the  soil  can  be  studied  by  consider'ng  an 
energy  balance  fur  that  material.  The  equation  expressing  the  balance  of 
mechanical  and  thermodynamic  energy  for  the  motion  of  a  general  continuum  can  be 
written  as  (Truesdell  and  Toupin  1960) 


A(KF)  =  j  Pdt  +  j  Qdt  -  A( IE) 
0  0 


(3.5) 


with 


P  =  (j)  t  *  v  da  +  IpB’vc 

Js  •'R 

Q  =  I  o  r  Jv  -  (j)  q  ’  n  da, 

^  -*5 


(3.6) 


(3.7) 


and  where  the  symbols  are  defined  as 


a(KE):  increment  of  kinetic  energy 

A(IE):  increment  of  internal  energy 

P:  rate  of  work  due  to  forces 

Q:  rate  of  he?' ing 

t:  surface  actions 

v:  velocity  of  material  particles 

p:  mass  density 

b:  body  force  per  unit  mass 

r:  internal  heating  rate  per  unit  mass 

q:  heat  flux  vector 

da:  area  differential  tlement 
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dv : 

volume  differential 

element 

dt: 

time  differential 

+ 

n : 

unit  vector  normal 

to  da 

and  R  is  any  material  region  enclosed  by  the  closed  surface  S. 

In  the  present  application,  these  equations  are  applied  to  the  region 
of  soil  material  excavated  to  create  the  crater,  the  region  labeled  as  R  in  Fig. 
3.2.  The  coupling  of  the  source  energy  is  via  a  pressure  pulse  transmitted 
across  the  source-soil  and  the  air-soil  interface,  denoted  as  the  surface  in 
Fig.  3.2.  Energy  transport  via  radiation  or  heat  conduction  is  inconsequential 
for  high-explosive  sources:  consequently,  Q  -  0.  The  rate  of  mechanical  work 
across  this  interface  is  given  by 


P 


1 


(3.8) 


This  rate  of  work,  integrated  over  the  duration  of  the  process,  is  that  fraction 
of  the  total  energy  that  does  work  on  the  excavated  material.  For  near-surface 
events,  much  of  the  energy  is  released  to  the  atmosphere.  Other  fractions  of 
the  total  energy  lost  are  in  the  shock  compression  and  subsequent  adiabatic 
unloading  of  the  ground  material  and  in  the  internal  shearing  of  the  material 
being  retained  as  internal  energy.  Recognizing  these  loss  mechanisms,  it  is 
assumed  that 

t 

P1  dt  -  A(IF)  =  nE  (3.9) 
o 

where  E  is  the  total  energy  of  the  explosive  and  n  is  that  fraction  that 
contributes  to  the  crater  excavat'on.  Consequently,  if  denotes  the  interface 
between  the  excavated  and  the  remaining  soil  (the  crater  boundary  interface, 
Fig.  3.2),  and  work  at  the  soil-air  interface  is  ignored,  then 


dt 


(3.10) 


which  states  that  the  energy  nE  produces  an  increment  in  kinetic  energy  and  does 
work  against  the  body  forces  and  along  the  crater  interface.  The  net  change  in 
kinetic  energy  is  zero  at  the  termination  of  the  process.  The  work  against  the 
gravity- induced  body  forces  can  be  written  as  an  increment  of  potential  energy. 
Thus , 


nE  =  a {PE)  +  wy  (3.11) 


where 

MPE)  =  - 

/ 

Jo 

f  pS  • 

•'R 

is  the 

increment  of  potential 

energy, 

and 

My  *  ^ 

u 

\  *  v  da 

Jo 

n 

l  J 

is  the  work  done  on  the  excavated  material  at  the  crater  boundary.  This  balance 
equation  states  that  the  energy  nE  goes  into  two  parts:  the  change  in  potential 
energy  required  during  the  excavation  and  the  work  done  as  the  excavated 
material  shears  along  the  crater  boundary.  Since  this  deformation  requires 
overcoming  the  strength  of  the  soil  material,  it  is  referred  to  as  the  work  done 
against  the  material  strength.  This  fundamental  energy  balance,  but  without  the 
n  coupling  factor,  was  attributed  to  Charters  and  Summers  (1959)  and  by  Gault 
and  Wedekind  (1977). 

The  exact  value  of  the  two  terms  on  the  right  in  Eq.  3.11  wil1 2  depend 
on  the  details  of  the  flow  field  during  the  cratering  process  and  on  the  final 
crater  configuration.  However,  by  using  certain  simplifying  assumptions  and 
considering  the  dimensions  involved,  the  form  of  these  two  terms  can  be  deduced. 

In  particular,  it  is  now  assumed  that 


1)  the  crater  is  hemispherical  with  radius  r; 

2)  the  work  at  the  crater  interface  is  proportional  to  the  material 
strength  Y ; 


3)  the  material  strength  Y  is  given  by  the  Mohr-Coulomb  model  (Eq. 

3.1),  where  c  is  the  cohesion,  P  the  confining  pressure,  and  ♦  the 
angle  of  internal  friction;  and 

4)  A  measure  of  the  confining  pressure  P  is  the  lithostatic  pressure 
p:  at  the  bottom  of  the  crater. 

If  the  strength  Y  (Eq.  3.1)  is  substituted  for  the  stress  t  in  Eq. 
3.12,  two  terms  result.  The  first  is  proportional  to  the  cohesion  c  and  the 
second  to  tan  ♦.  An  integration  over  the  crater  surface  S2.  taken  here  to  be  a 
hemisphere  for  convenience,  will  give  the  form: 

Wy  =  C^cr3)  +  C2  (pgr4  tan  *) ,  (3.13) 

where  and  C2  are  two  numerical  constants  which  depend  upon  the  crater 

geometry.  The  form  of  Eq.  3.13  is  general  and  applies  to  any  axi-symmetric 

family  of  crater  shapes  where  the  depth  is  a  given  fraction  of  the  radius. 

The  first  term  results  from  the  cohesion  c  acting  on  a  surface  area 
? 

proportional  to  r  and  shearinq  along  a  distance  proportional  to  r.  The  second 

term  is  due  to  a  confining  pressure  proportional  to  Pgr  acting  on  the  same 

2 

surface  that  is  proportional  to  r  and  sheared  over  a  distance  proportional  to 
r. 

Likewise,  the  potential  energy  term  is  given  by 

A(PE)  =  C3  pgr4  ,  (3.14) 

3 

where  the  excavated  volume  is  proportional  to  r  *  the  elevation  change  is 
proportional  to  r  and  the  constant  C3  is  determined  by  the  geometry  of  the 
crater  shape.  Consequently,  the  energy  balance  has  the  general  form 

nE  =  C-|  (cr3)  +  C2  (pgr4  tan  $)  +  C3(pgr4).  (3.15) 

The  values  of  ,  C2»  and  C3  also  depend  on  the  flow  field  leading  to  the  crater 
excavation.  For  example,  Gault  and  Wedekind  (1977)  give  a  similar  form,  but 
without  the  n  factor,  or  the  tan  4>  term,  and  with  values  for  C-|  and  C3  based 
upon  additional  simplifying  assumptions. 

The  energy  of  the  explosive  given  in  terms  of  the  charge  radius,  a,  is 
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(3.16) 


E  =  (4«/3)  a3  «Qe. 

The  numerical  factor,  4«/3,  can  be  combined  with  the  other  coefficients  in  Eq. 
3.15  when  Eq.  3.16  is  used  to  get 

n  ’  ^  ^  ^  +  ^8^  ^  '“2  tan  *  +  C3).  (3.17) 


The  dimensionless  groups  defined  in  Eqs.  3.2  can  be  rewritten  as 

(3.i 8a) 
(3.18b) 
(3.18c) 


*  =  <L 
4  °Qe 

■  tan  4. 

Therefore  Eq.  3.17  can  be  written  as 


(3.1 8d ) 
(3.1 8e ) 


n  =  (»4)  (l»y)  +  (*2)  {*y)  '  ("3)' 


+  K  o ) . 


(3.19) 


The  factor  *g  =  °/6.  which  occurs  to  the  negative  1/3  power,  does  not  vary 
significantly  over  the  ranqe  of  experiments  considered  here.  As  a 
simplification,  it  is  assumed  to  be  constant  and  is  henceforth  included  in  the 
K,p  and  Kg  coefficients  giving 

n  =  K]  wv  +  (K2  +  K3)  *2  "y4/3.  (3.20) 

Equation  3.20  can  be  used  to  express  wy  as  a  function  of  the  other 
parameters  *g  and  ir^,  recognizing  that  the  coupling  factor  n  is  expected  to 
depend  upon  "2*  "3  an(^  "4'  Only  after  specifying  this  dependence  is  it  possible 
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to  determine  the  explicit  form  of  the  solution  to  this  equation.  However,  it  is 
instructive  to  consider  certain  limiting  cases.  For  example,  consider  a  series 
of  experiments  in  a  medium  having  zero  cohesion  c.  In  this  case,  ^  =  0,  is 
a  constant  and  Eq.  3.20  reduces  to 


(3.21) 


For  a  series  of  experiments  in  a  cohesionless  material  for  a  given  explosive, 
and  assuming  that  the  efficiency  n  is  independent  of  *2, 


V 


e3/4  q-3/V 


(3.22) 


In  this  case,  at  fixed  gravity,  each  crater  linear  dimension  varies  with  the 
one-quarte^  root  of  the  energy,  often  referred  to  as  quarter-root  scaling. 
However,  experimental  results  (Piekutowski ,  1975;  Schmidt  and  Holsapple,  1978a, 
1980)  indicate  that  the  coupling  efficiency  n  depends  upon  Qp  and  g;  hence,  n  is 
not  independent  of  i^*  Furthermore,  it  can  also  be  expected  to  depend  on  "3* 
Consequently,  quarter-root  scaling  has  not  been  observed  for  half-buried  high 
explosive  charges,  even  in  cohesionless  materials;  whereas  for  hypervelocity 
impact  into  water,  Gault  (1978)  obtained  cratering  data  that  approaches  quarter- 
root  seal ing--l /3.83. 

A  second  special  case  arises  when  the  term  with  *2  is  small  compared  to 
the  first  term  in  Eq.  3.2..  This  occurs  if  either  g  or  the  charge  energy  E  is 
sufficiently  small,  and  t  material  has  non-zero  cohesion.  Then,  noting  that 
1:4  is  constant  for  a  giver  soil-explosive  combination,  the  cratering  efficiency 
is  proportional  to  the  coup. inq  factor, 


IT  a 

V 


n . 


(3.23) 


Furthermore ,  if  n  is  independent  of  «2  and  then  the  cratering  efficiency  is 
constant  ana 

V  «  E  ,  (3.24) 


resulting  in  the  cube-root  scaling. 


In  the  general  case  for  a  soil  with  finite  cohesion,  allowing  the 
coupling  efficiency  n  to  be  a  variable,  the  solution  for  the  cratering 
efficiency  will  be  determined  by  all  the  terms  in  Eq.  3.20.  That  is 

%  =  G  [n,  *4,  (*5  +  k 1 ) J  (3.25a) 

where  k]  =  K3/K2-  (3.25b) 

A  study  of  the  experimental  cratering  data  for  a  given  soil-explosive 
combination  shows  that  the  results  follow  a  consistent  trend  with  increasing 
charge  weight.  At  small  values  of  w 2  (small  charge  size),  the  cratering 
efficiency  i»v  approaches  a  constant  (horizontal  asymptote).  As  tt2  increases,  *v 
decreases,  approaching  a  different  asymptote  for  large  What  is  desired  is  a 
functional  relationship  that  fits  this  general  trend.  In  Eq.  3.25a,  the 
variables  *2  and  tr g  occur  in  a  particular  combination.  Consequently  the 
cratering  efficiency  can  be  expected  to  depend  on  this  combined  parameter  and 
separately  on  ^4.  This  additional  parameter  which  is  in  principle  a 

completely  distinct  variable,  is  included  as  a  simple  sum 

w2  =  w2  (ff5  *  ^1^  *  ^2  "4  (3.26a) 

=  %  (-f)1/3  Ctan  *  +  +  k2  .  (3.26b) 

we  MMe 

This  single  variable  w2  includes  the  strenqth  measures  c  and  tan  0,  as  well  as 
the  gravity-seal ed-yi el d  parameter  "2  and  the  specific  energy  measure  of  the 
explosive  Q  Whether  it  alone  can  be  used  to  account  for  dependences  on  *2, 
and  separately  can  only  be  determined  from  experiments.  Assuming  that  it 
can,  then 


"v  =  H  (?2)  (3.27) 

for  all  soils  and  explosives. 

For  the  specific  case  of  the  dry-Ottawa-sand  experiments,  the  cohesion 
c  is  zero.  Therefore,  *4  is  zero  and  Eq.  3.26a  becomes 

*2  =  (tan  0  +  kj  )t>2  (3.28) 
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where  the  sum  (tan  4>  +  k-j)  is  a  constant.  The  variable  is,  therefore,  only  a 
constant  multiple  different  from  the  gravity-scaled-yield  parameter  *2.  As  shown 
in  Fig.  3.1,  the  sand  data  are  well  represented  by  a  power  law  of  the  form 
(Schmidt  and  Holsapple  1978a,  1980) 

if v  =  0.194  tt2'0-472  .  (3.29) 
In  terms  of  the  i?2  variable,  Eq.  3.29  can  also  be  written 


7T 

V 


0.194 


tan  ♦  + 


-0.472 


(3.30) 


The  measured  value  of  $  for  this  material  was  35°  giving  tan  4>  =  0.70. 
Experimental  results  for  materials  with  different  values  for  tan  determine  the 
empirical  constant  k^.  A  value  of  0.1  was  chosen  to  be  consistent  with  the 
various  soils  tested.  The  factor  k-,  is  a  measure  of  the  relative  importance  of 
the  potential  energy  term  compared  to  the  Mohr-Coulomb  confining-pressure 
strength  term  in  the  energy  balance  (Eq.  3.15).  The  fact  that  k-|  is  relatively 
small  compared  to  tan  4>  for  granular  materials  indicates  that  the  dominant 
influence  of  increased  gravity  or  increased  charge  size  in  crater  formation  is 
not  the  increasing  potential  energy  contribution,  but  the  increased  lithostatic 
pressure.  Only  for  a  material  with  zero  angle  of  internal  friction  will  this 
dependence  on  increasing  size  be  dominated  by  the  potential  energy  contribution. 

For  k^  =  0,1,  Eq.  3.30  becomes 


n 

v 


0.194 


(0.7  +  0.1) 


-0.472 


-0.472 


(3.31) 


giving  the  form  of  the  function  H  in  Eq.  3.27 

(*2)  =  0.174  ii2'0-472  (3.32) 

The  usefulness  of  the  particular  combination  of  terms  defining  *2, 
(Eqs.  3,26)  is  shown  in  Fig.  3.3  for  various  soil-explosive  pairs.  For  all 
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materials  tested,  »  versus  the  parameter  can  be  fitted  by  a  single  straight 
line  on  a  log- log  plot. 

The  evaluation  of  the  ^  parameter  requires  values  for  the  strength 
properties  of  the  various  soils,  as  well  as  the  material  properties  of  the 
explosives.  Experimentally  determined  values  for  the  cohesion  c  ana  the  anqle 
of  internal  friction  <l>  were  measured  using  static  triaxial  tests  (see  Appendix 
B).  These  values  are  listed  in  Table  3.1  for  dry  Ottawa  sand,  KAFB  desert 
alluvium,  and  the  oil-base  clay.  Table  3.2  lists  the  values  of  specific  energy 
Qg  and  mass  density  5  for  the  various  charges  used  in  this  study.  The  three 
sizes  of  PETN  charges  have  different  Qe,  because  a  fixed  amount  of  silver  azide 
initiator  was  used  in  each  of  the  various  sizes.  The  value  of  the  specific 
energy  Qg  for  each  composite  charge  was  calculated  by  dividing  the  energy 
release  due  to  both  the  silver  azide  and  the  PETN  by  the  sum  of  their  masses. 

A  value  of  unity  for  the  other  empirical  constant  in  Eq.  3.26,  was 
determined  by  trial  and  error.  The  final  results  are  shown  as  a  curve  of 
versus  for  the  three  materials  and  the  two  different  explosives  in  Fig.  3.3. 
In  contrast  to  Fig.  3.1,  the  combined  gravity-yield-strength  parameter  ^ 
accounts  for  the  differences  in  gravity,  size,  strength  and  charge  properties, 
with  particularly  good  agreement  for  the  sand  and  alluvium  experiments.  All  the 
data  points  lie  near  a  single  straight  line.  Only  the  clay  data  shows 
appreciable  scatter,  which  is  discussed  below.  These  results  provide 
experimental  justification  of  the  particular  combination  of  parameters  defining 
the  *2  parameter. 

It  is  instructive  to  rewrite  the  function  ity  =  H  ("2)  in  terms  of  the 
gravity-scaled-yield  parameter  *,,,  and  the  strength  parameters  and 
Substituting  Eq.  3.26a  into  Eq.  3.32  gives 

l/n  r  n  w  l/3  -0.472 

^  =  je.  =  0.174  C-fg-  "  §  (?)  (tan  *  +  0.1)]  .  (3.33) 


Curves  generated  from  this  equation  for  each  different  combination  of 
soil  cohesion  and  charge  properties  are  shown  in  Fig.  3.4.  It  is  seen  that  this 
equation  does  reproduce  the  observed  experimental  trends  and  values.  The 
agreement  is  particularly  good  for  the  sand  and  the  alluvium  data  and  correctly 
matches  the  shape  for  both  explosive  types  in  the  clay  which  exhibits  some 
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Table  3.1.  Soil  Properties. 


Cohesion 

Cratering  Medium 

p 

(dyne/cm  ) 

Dry  Ottawa  Sand 

0.0 

Dry  Desert  Al 1 uvium 

1.7  x  105 

Model inq  Cl  ay 

1.1  x  105 

Interna!  Friction  Angle  Density 


(degree) 

3 

(gm/cm  ) 

35 

1.80 

32 

1 .60 

1.2 

1.53 

Table  3.2.  Explosive  Properties 


Explosive 

Type 

Mass 

(<P) 

Energy 

(erg) 

Specific  Energy 
(erg/gm) 

Density 

3 

( gm/cm  ) 

PETN 

0.49 

2.24  x  1010 

4.54  x  1010 

1.93 

PETN 

1.34 

7.18  x  1010 

5.34  x  1010 

1.78 

PETN 

4.08 

2.27  x  1010 

5.56  x  1010 

1.73 

PbNg 

1.70 

2.24  x  1010 

1.32  x  1010 

3.10 
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scatter.  Possible  causes  of  this  scatter  include  rate  effects  which  violate  the 
similarity  requirements  at  fixed  *£»  sample  preparation,  aging  of  the  clay  and 
dependence  of  material  properties  upon  temperature  which  varied  slightly  from 
shot  to  shot  (see  Table  2.3).  Additional  experiments  are  required  to  better 
define  the  clay  behavior  leading  to  the  observed  scatter.  Nevertheless  the  mean 
values  are  in  good  agreement. 

For  small  values  of  *2  (sma^  charge  mass  W),  *v  aporoaches  a  constant 
value  given  by 


r  -0.472 

%  =  0J74  (3-34) 

and  is  therefore  determined  by  the  cohesion  and  density  of  the  cratering  medium 
and  the  specific  energy  of  the  explosive.  Equation  3.34  agrees  with  the 
differences  in  the  observed  asymptotes  for  the  various  combinations  tested. 
Only  if  c  is  zero  will  this  cube-root  regime  not  be  observed,  as  is  the  case  for 
the  dry  Ottawa  sand  data. 

For  sufficiently  large  values  of  charge  size,  the  * ^  term  in  Eq.  3.26a 
will  dominate.  In  this  large-charge-size  regime,  the  curves  generated  from  Eq. 
3.33  are  asymptotic  to  the  equation 

n  u  V3  -0.472 

=  0.174  [^-  (|)  (tan  *  +  0.1)]  (3.35) 

which  is  independent  of  the  cohesion  of  the  material.  This  is  a  straight  line 
on  a  log-log  plot  of  ir  versus 

3-3  GENERALIZATION  TO  OTHER  SOIL  TYPES 

Figure  3.5  shows  the  expected  behavior  for  a  generic  material  having 
both  finite  cohesion  and  finite  angle  of  internal  friction.  The  horizontal 
asymptote  described  above  is  shown,  as  well  as  the  asymptote  for  large  values  of 
charge  size.  The  transition  value  of  corresponding  to  the  intersection  of 
these  two  straight  line  asymptotes,  is 


c 

*2  =  p Q  (0.1  +  tan  ♦)  ’ 


(3.36) 


37 


Fig.  3.5  General  form  of  cratering  efficiency  versus  gravity-scaled  yield  for  a  typical  soil. 


which  is  determined  by  the  cohesion,  the  angle  of  internal  friction,  the  soil 
density  and  the  specific  energy. 

Using  the  definition  of  *2»  the  critical  charge  mass  corresponding  to 
this  transition  is 


W  6  [  pg  (tan  <t>  +  0. 1  )J  *  (3.37) 

A  curve  depicting  the  transitional  charge  size  for  TNT  in  terms  of  the 
strength  parameter  c/p  (tan  $  +  0.1)  for  various  values  of  gravity  g  is  shown  in 
Fig.  3.6.  The  strength  ranges  for  common  soils  are  indicated  on  this  figure  in 
a  very  approximate  way.  This  figure  indicates,  for  example,  that  at  1  G  the 
transitional  charge  mass  for  alluvial  materials  is  in  the  range  of  a  ton  to  a 
kiloton  of  TNT. 

These  results  explain  some  of  the  apparent  discrepancies  observed  in 
explosive  cratering  tests.  In  alluvium,  for  example,  small-scale  tests  do  not 
extrapolate  to  large-scale  field  tests,  as  can  be  seen  based  upon  the  strength 
model  given  here.  Small  tests  in  alluvium  at  1  G  are  in  the  cohesion-dominated, 
cube-root  regime.  Existing  large-scale  field  tests  are  around  the  transition 
charge  size.  These  field  data  (Chabai ,  1965)  are  compared  with  the  calculated 
behavior  for  TNT  explosive  (Q  =  4.19  x  1010  erg/gm)  in  desert  alluvium  as  shown 
in  Fig.  3.7.  The  predicted  crater  volume  is  based  entirely  upon  the  centrifuge 
experimental  results,  using  charge  sizes  on  the  order  of  a  few  grams  of  either 
PETN  or  PbNg.  The  5000-lb  shot  lies  right  on  the  predicted  curve.  There  is 
some  scatter  amonq  the  three  256-lb  shots,  but  the  mean  value  lies  on  the  curve. 

Another  interesting  result  is  a  comparison  of  the  clay  and  the  sand 
cratering  efficiencies.  At  small  n  values  (small  charge  size)  the  clay  volumes 
were  substantially  below  those  for  sand.  This  is  consistent  with  the  idea  that 
sand  is  "strengthl ess"  and  clay  has  finite  cohesion.  For  large  (large  charge 
size)  the  cratering  efficiency  in  clay  is  greater  than  in  dry  sand.  In  this 
regime,  the  cratering  efficiency  is  dominated  by  the  angle  of  internal  friction, 
which  introouces  an  effective  strength  due  to  the  large  lithostatic  pressure. 
Clay  has  a  near-zero  angle  of  internal  friction.  Consequently,  for  large 
confining  pressure,  the  dry  sand  is  "stronger"  than  the  clay,  which  explains  the 
observed  crater  volumes.  This  general  behavior  is  in  qualitative  agreement  with 
calculated  results  of  O'Keefe  and  Ahrens  (1578)  who  used  a  Mohr-Coulomb  strength 
model  with  a  high-pressure  cutoff. 
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SECTION  4 

CRATER  RADIUS  AND  CRATER  DEPTH 

Significant  variations  in  crater  shape  versus  scaled  energy  were 
observed  for  the  different  soil-explosive  combinations  tested.  These  variations 
include  differences  in  overall  shape,  as  measured  by  the  aspect  ratio  (ratio  of 
the  radius  to  depth),  and  more  detailed  features,  such  as  terraces  and  central 
mounds.  For  the  purposes  of  this  report  only  the  differences  in  radius  and 
maximum  depth,  are  considered. 

Nondiinensional  forms  for  the  crater  radius  r  and  the  crater  depth  h  are 

1/3 

(4J) 

1/3 

=  h  (-£)  (4.2) 

All  the  results  obtained  to  date  for  nr  and  for  wh  for  dry  sand,  clay, 
and  4-percent  moisture  alluvium  at  zero  depth  of  burst  are  shown  in  Figs.  4.1 
and  4.2  as  a  function  of  the  gravity-scaled-yield  parameter  *£•  The  data  points 
for  dry  Ottawa  sand  are,  as  given  previously  (Schmidt  and  Holsapple  1978a), 
along  the  curves  described  by 

=  0.766  *2‘0'159  (4.3) 

and 

irh  =  0.164  n2"°*164  (4.4) 

which  give  very  good  straight-line  fits  using  these  logarithmic  scales. 

For  the  alluvium  data  points,  considering  only  those  with 
(approximately)  4-percent,  moisture,  the  acquisition  of  additional  data  has  led 
to  an  interpretation,  incorporating  the  cohesion  of  the  medium,  different  than 
that  previously  qiven  (Schmidt  and  Holsapple  1978a).  The  data  for  the  radius 
has  the  same  qualitative  trend  as  does  the  crater  volume  data.  A  functional  form 
for  a  quantitative  model  is  qiven  below.  The  trend  of  the  depth  data,  however, 
is  not  even  qualitatively  the  same  as  that  for  either  the  volume  dr  the  radius. 
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Fig.  4.2  Depth  parameter  versus  gravity-scaled  yield  for  various  soil -explosive  combinations. 


An  attempt  was  made  to  model  both  the  radius  and  depth  data  using  the 
same  strength  model  developed  for  the  volume  data  (Section  3).  However,  the 
data  do  not  appear  to  be  consistent  with  this  interpretation.  Some  reasons  for 
this  are  discussed. 

The  strength  model  for  volume  was  based  on  the  energy  available  and 
upon  the  work  necessary  to  excavate  the  crater.  However,  the  dependence  of  the 
work  required  on  the  specific  shape  of  the  crater  was  not  considered.  Thus,  the 
theory  is  expected  to  be  more  suitable  for  the  volume,  than  for  shape  variables 
such  as  radius  and  depth.  Since  the  volume  is  approximately  proportional  to  the 
radius  squared  times  the  depth,  the  radius  is  the  dominant  dimension  determining 
the  volume.  Consequently,  the  radius  behavior  might  be  expected  to  exhibit  the 
same  trend  as  the  volume  and,  perhaps,  the  depth  to  a  lesser  degree. 

It  is  clear  that  slope  stability  considerations  may  rule  out  certain 
crater  shapes  as  either  the  explosive  size  or  gravity  is  increased.  For 
example,  McKinnon  and  Melosh  (1978)  have  shown  for  craters  with  vertical  walls 
and  flat  bottoms  that  there  is  a  maximum  depth,  due  to  the  initiation  of 
slumpinq.  The  value  depends  upon  the  parameters  ogh/c  and  tan  4).  Subsequent 
work  by  Melosh  and  McKinnon  (1979)  indicates  that  craters  of  parabolic  shape 
also  cannot  exceed  some  maximum  depth, 
observations  of  James  (1977,  19/8). 

From  the  following  equality 

^  =  h(§)1/3  a-  (H-)1/3 

_  / /  1 /3  v 

"  ^ ^ 4  *3  % 

a  limiting  value  on  ogh/c  implies  a  limit  on  the  crater  depth  variable  as  a 
function  of  for  a  given  soi 1 -expl osi ve  combination  which  determines  tt4  and 

*3* 


This  trend  is  consistent  with 
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e  ,fiJ/3 


n 


c  p 


(4.5a) 

(4.5b) 


"h  1 


<  K 


1/3  -1 

w4  w3 


(4.6) 


The  constant  K  depends  on  tan  i>  and  crater  shape  profile.  If,  to  first 
order,  it  can  be  assumed  that  slumpinq  due  to  stability  limits  does  not  change 
the  crater  volume,  the  prediction  of  crater  volume  only  requires  a  consideration 
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of  the  energy  balance.  The  radius  is  influenced  to  some  degree  by  stability 
considerations,  and  the  depth  and  the  aspect  ratio  (r/h)  may  be  dominated  by 
stabil ity. 

The  centrifuge  does  give  a  tool  for  direct  simulation  of  large-scale 
craters.  For  the  soils  and  explosives  over  the  range  of  scaled  energy  tested, 
the  data  points  in  Figs.  4.1  and  4.2  give  predictive  values  for  both  radius  and 
depth.  For  granular  materials,  like  sand  and  alluvium,  wr  versus  ^  can  be 
explained  with  a  functional  form  like  that  obtained  for  the  volume.  The  clay 
craters  were  significantly  different  in  shape,  and  the  theory  generated  from 
energy  balance  considerations  alone  was  not  sufficient  to  explain  the  data. 

At  small  values  of  the  gravity-scaled  size  parameter  the 

variation  for  the  alluvium  is  essentially  constant  (independent  of  This  is 

a  regime  of  cube-root  scaling  where  all  the  crater  dimensions  are  determined  by 

the  cohesion  of  the  material.  The  dry  sand,  which  has  essentially  no  cohesion, 

does  not  show  this  behavior. 

As  the  variable  it?  is  increased,  the  results  become  more  influenced  by 
lithostatic  pressure.  As  discussed  in  Section  3,  this  is  due  to  two  factors: 
the  increased  strength,  due  to  the  lithostatic  pressure  and  the  angle  of 

internal  friction;  and  the  additional  work  that  must  be  done  to  excavate  the 
crater  as  the  size  or  the  gravity  increases.  For  the  volume  variable  nv,  the 
data  show  a  transition  regime,  curving  from  the  constant  it  asymptote  for  small 
it.?  values  to  a  sloped  straight-line  asymptote  (on  log-log  plot)  for  large 
values.  The  radius  parameter  n  for  alluvium  shows  the  same  general  trend,  but 
the  exact  form  of  the  function  that  fits  the  data  is  different. 

As  given  above,  the  effects  of  strength  for  the  various  materials 
considered  can  be  measured  by  a  general ization  of  the  combined  strength- 
gravity-scaled  yield  parameter  ii^  given  by  Eq.  3.26b  as  follows 

*2  =  ^A1  tan  *  +  A2^  +  A3  pQ  ‘  (4*7^ 

There  are  three  terms  in  this  parameter.  The  A^  tan  *  term  is  a  measure  of  the 
work  against  the  material  strenqth  increase  due  to  confining  pressure.  The  A^ 
term  is  a  measure  of  the  work  against  gravity,  and  the  A^  term  measures  the  work 
to  overcome  the  cohesion  during  the  crater  excavation.  The  valiies  of  each  of 
the  constants  determine  the  relative  importance  of  each  of  these  terms. 
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For  the  volume  data,  was  taken  as  unity,  A2  was  taken  to  be  0.1,  and 
A j  as  1.0,  as  discussed  in  Section  3.  However,  for  the  radius  data,  different 
values  are  seen  to  be  more  appropriate.  This  is  not  surprising,  since  the 
different  work  terms  vary  in  importance  in  determining  radius,  as  compared  to 
the  depth  or  the  volume. 

As  shown  in  Fig.  4.1,  the  values  for  *  for  alluvium  and  dry  Ottawa 
sand  all  converge  at  large  values  of  it2.  This  indicates  that  there  is  little  or 
no  dependence  on  tan  $,  since  this  strength  contribution  determines  any  varia¬ 
tion  in  dependence  on  *2  for  large  values  of  «2.  Thus,  for  the  radius  data,  A^ 
was  taken  to  be  zero  and  ii2  can  divided  by  A2  giving 


«  i  (V  +  A 

"2  QeV  A4  pQe* 

where  now 

A 4  *  A3/A2  . 


(4.8a) 

(4.8b) 


To  find  the  best  value  of  A^,  the  cohesion-determined  asymptotes  at 
small  values  of  *2  can  be  used.  The  functional  form  expected,  based  on  the 
results  of  Section  3,  is 
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(4.12b) 


a  =  0.159. 

For  materials  other  than  sand  this  result  can  be  generalized  by 
substituting  the  complete  definition  of  ?  (Eq.  4.8a)  into  Eq.  4.11,  giving 


a  W  1 /0  r  -0.159 

v  =  0.766  [^-  (|)I/J  +  A4  ^-]  .  (4.13) 

The  constant  A4  is  determined  by  the  value  of  the  cohesion-dominated  asymptotes 
at  small  values  of  Setting  ^  =  0, 


0.766  (A4r0J59  (^-)"°'159- 


(4.14) 


Alluvium  data  for  tr  at  small  values  are  taken  from  Piekutowski  (unpublished 
data)--shot  numbers  UDRI-646  through  UORI-651  in  Table  2.2.  For  these  shots,  A^ 
=  0.5  gives  a  good  fit.  This  completely  fixes  the  dependence  of  „  on  the  other 
variables  and  the  final  result  is 


n  u  1/3  .  -0.159 

V  =  0.766  [§  (f)  +  0.5^1  (4.15) 

6  pUp 

The  curves  generated  from  this  equation,  together  with  the  data  points  for  dry 
Ottawa  sand  and  alluvium,  are  shown  in  Fig.  4.3. 

As  a  final  measure  of  crater  shape,  the  aspect  ratio,  r/h,  will  be 
discussed.  As  suggested  above,  slope  stability  considerations  are  expected  to 
play  an  important  role  in  determining  values  of  this  crater  dimension.  The 
slope  stability  depends  upon  the  material  strength,  which  has  two  components: 
the  cohesion  and  the  angle  of  internal  friction.  For  a  material  with  zero 
cohesion  and  non-zero  angle  of  internal  friction  (such  as  the  dry  Ottawa  sand), 
stability  is  independent  of  gravity  and  size.  That  is,  stability  of  a  certain 
shape  at  small  scale  or  at  small  gravity  implies  stability  at  large  scale  or  at 
large  gravity.  This  follows  from  the  fact  that,  while  increased  size  or  gravity 
requires  greater  strength,  the  strenqth  increases  as  the  size  or  gravity 
increases  due  to  the  increased  lit hos t a t i c  confining  pressure  (i.e.,  strength 
proportional  to  ogh  tan  <j>). 
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strength  model  for  radius  parameter  with  centrifuge  experimental  results 


On  the  other  hand,  the  stability  of  a  material  with  small  angle  of 
internal  friction  but  large  cohesion  will  depend  both  on  the  scaled  size  and  on 
gravity.  Consequently,  a  variation  of  shape  is  ant.'ipated  with  increasing 
si  ze . 

The  parameter  * ^  correlates  all  crater-size  data  very  well.  It 
accounts  for  differences  of  size,  gravity  and  source  energy  density.  However, 
for  shape  measures  such  as  the  aspect  ratio,  *2  may  not  be  the  most  appropriate 
independent  parameter.  For  example,  energy-density  effects  may  not  play  a  role 
in  crater  shape.  The  stability  arguments  given  above  suggest  that  for  a 
cohesive  material  the  dimensionless  parameter  pgh/c  may  be  more  significant. 
(It  should  be  noted  that  this  pi-group  is  not  independent  but  merely  a 
combination  of  the  original  set  as  shown  in  Eq.  4.5b.) 

Figures  4.4,  4.5  and  4.6  show  aspect  ratio  results  for  the  various 
soils.  Figure  4.4  is  a  composite  plot  showing  the  aspect  ratio  for  alt  three 
materials  versus  Figure  4.5  shows  the  alluvium  aspect  ratio  versus  pgh/c, 
and  Fig.  4.6  shows  the  same  for  clay.  The  expectations  are  supported  by  the 
data.  For  dry  Ottawa  sand,  the  aspect  ratio  is  essentially  independent  of 
scaled  size,  as  indicated  on  Fig.  4.4.  The  value  is  about  4.6  with  +.10  percent 
variation. 

For  the  alluvium,  Fig.  4.5  is  thought  to  be  more  meaningful  than  Fig. 
4.4.  In  Fig.  4.5,  with  the  stability  parameter  pgh/c  as  the  abscissa,  the 
aspect  ratio  shows  as  a  gradual  decrease  (more  hemispherical  shape)  as  the 
parameter  increases.  For  comparison,  the  equivalent  1-G  crater  depth  is  shown 
on  tho  top  as  a  second  abscissa  scale.  This  trend  is  not  expected  to  continue 
to  ever- i ncreasi ng  crater  sizes.  For  very  large  sizes,  slope  stability  dictates 
that  the  crater  walls  can  never  qet  steeper  than  the  angle  of  friction  for  the 
material.  A  lower  limit  on  r/h  should  be  given  by  the  reciprocal  of  tan*.  For 
alluvium,  this  limit  value  is  r/h  =  1  .6  as  shown  in  Fig.  4.5.  With  increasing 
crater  depth  the  aspect  ratio  is  expected  to  become  asymptotic  to  this  limit. 
Interestingly,  the  TEAPOT  ESS  crater,  whose  depth  was  20.4  in,  had  the  smallest 
aspect  ratio  of  all  the  large  NTS  alluvium  craters,  a  value  equal  to  1.62. 

Figure  4.6  shows  the  clay  results  as  a  plot  of  aspect  ratio  versus  the 
stability  parameter.  For  small  values  of  pgh/c,  the  craters  are  nearly 
hemispherical .  As  the  stability  parameter  increases,  there  is  a  gradual 
increase  in  the  aspect  ratio.  Only  one  data  point,  shot  25- X  has  a  uniquely 
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il-base  clay. 


greater  aspect  ratio  (flatter  crater) .  It  also  has  the  largest  value  of  the 
stability  parameter,  which  probably  accounts  for  its  unique  flat  profile. 
McKinnon  and  Melosh  (1978)  have  found,  for  "transient"  craters  with  a  cake-pan 
shape,  that  a  stability  limit  occurs  at  a  definite  value  of  pgh/c  for  a  given 
material  with  a  given  tan  For  cake-pan-shaped  craters,  they  found  the  value 
to  be  about  6.0.  For  the  present  experiments,  craters  seem  to  retain  a  lower 
value  of  aspect  ratio,  up  to  a  value  of  pgh/c  equal  to  approximately  18.  This 
is  equivalent  to  a  value  of  pgh  equal  to  2.1  x  10  dyne'rm1 .  shot  25- X  is  the 
only  shot  over  this  threshold. 

A  check  on  the  existence  of  a  stability  limit  based  on  Dgh/c  was  made 

as  an  adjunct  to  shot  25-0.  This  crater  was  formed  at  low  G  with  a  resulting 

pgh/c  equal  to  0.8.  Subsequent  to  the  cratering  event,  the  qravity  was 

increased  to  518  G  for  100  seconds,  which  woi  J  have  resulted  in  a  pgh/c  equal 

to  42  had  the  crater  shape  remained  constant.  However,  the  crater  slumped 

considerably  as  can  be  seen  in  Figure  4.7.  The  final  aspect  ratio  was  2.1, 

about  the  same  as  that  for  crater  25-X  which  was  formed  at  high  G.  The  final 

slumped  values  of  pgh/c  for  25-0  was  equal  to  23,  and  Dgh  was  2.5  x  10 
2 

dynes/cm  .  This  gives  compelling  evidence  that,  for  this  clay,  all  large  volume 

p 

craters  will  have  a  limiting  value  of  pgh  equal  to  2.5  x  10  gm/cm-sec  .  At  1 

3 

G,  with  p  =  1,53  gm/cm  ,  the  limiting  depth  is  about  16  meters  for  all  charge 
si zes . 

Crater  25-X,  formed  at  hiqh  G,  was  also  subjected  to  the  518-G 
environment  for  100  seconds  to  provide  a  control  point.  As  can  be  seen  in  Fig. 
4.7,  it  underwent  negligible  shape  chanqe  suggesting  that  a  stable  shape  was 
formed  upon  excavation.  This  was  confirmed  when  both  craters,  25-0  and  25-X, 
were  cut  in  two  using  a  thin  wire  to  reveal  the  crater  cross  sections. 

The  low-G  crater,  25-0,  showed  marked  evidence  of  flow  and  collapse  of 
the  transient  crater  formed  at  low  G,  whereas  the  high-G  crater,  25-X,  showed  no 
evidence  of  a  transient  crater  or  any  significant  flow  or  slumping.  Figure  4.8 
is  a  comparison  of  the  two  cross  sections.  They  show  the  same  approximate  final 
shapes,  but  were  arrived  at  by  very  different  processes. 

In  this  experiment,  using  the  oil-base  clay  material,  the  crater  formed 
at  hiqh  G,  corresponding  to  a  1-G  field  event  of  564  metric  tons  of  PETN,  did 
not.  go  through  a  transient  stage  and  collapse. 
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Fig.  4.7  Slumping  due  to  subsequent  high-G  loading. 

(Also  see  Figs.  A23-A26,  Fig.  4.8  and  text. 


rig.  4.8. 


Comparison  or  crater  cross  sections.  Shot  25-0  formed  at 
10  G  shown  to  the  right.  Shot  25-X  formed  at  51 7  C.  shown 
to  the  left.  Both  craters  were  ^uhsequentl v  spun  at  518  G 
for  100  seconds.  Note  flow  pattern  and  structure  of  crater 
floor  in  10-G  crater.  (Also  see  riys.  A23-A26,  Tig.  4.7  and 
text. ) 
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SECTION  5 
MOISTURE  CONTENT 


Several  shots  were  designed  to  determine  the  effect  of  moisture  content 
on  cratering  phenomena  in  desert  alluvium.  Four  shots  with  moisture  contents 
from  less  than  0.7  percent  to  7.5  percent  were  made  in  addition  to  the  series  of 
shots  at  approximately  4-percent  moisture. 

Significant  differences  due  to  moisture  content  were  found.  Static 
tri axial  tests  were  performed  on  specimens  with  various  moisture  contents.  The 
strength  properties  were  found  to  depend  strongly  upon  moisture  content. 
Consequently,  it  is  expected  that  differences  in  cratering  behavior  should  be 
observed.  The  results  of  the  static  material-property  tests  are  summarized  in 
Table  5.1.  Two  values  of  cohesion  are  shown  for  each  case.  One  of  these  is 
based  on  the  zero-conf i ni ng-presoure  intercept  from  triaxial  tests.  The  other 
is  obtained  by  a  direct  shear  test. 

All  variable  moisture  tests  were  performed  at  a  gravity-scaled-yield 
parameter  " ^  about  equal  to  7  x  10"^.  This  corresponds  to  about  70  metric  tons 
of  TNT  at  1  G.  The  variation  of  cratering  efficiency  *  with  moisture  content 
at  this  *2  value  is  shown  in  Fig.  5.1  (also  see  Fig.  3.4).  The  dependence  of 
aspect  ratio  r/h  on  moisture  is  shown  in  Fig.  5.2. 

Although  an  explicit  correlation  to  measured  strengths  was  not 
attempted  since  all  the  data  were  for  a  fixed  scaled  size,  certain  qualitative 
behavior  is  noteworthy.  Both  the  crater  size  and  shape  depend  markedly  on 
moisture  content  at  this  scaled  size.  The  cratering  efficiency  decreases  as  the 
moisture  content  decreases,  with  the  largest  changes  occurring  as  the  alluvium 
becomes  more  dry.  This  is  probably  due  to  the  significantly  larger  angle  of 
internal  friction  noted  for  the  dry  alluvium.  At  this  scaled  size,  the  value  of 
^2  is  in  the  reqime  where  that  friction  angle  is  the  dominant  strength 
contribution. 

Since  all  samples  with  moisture  content  greater  than  about  2  percent 
have  about  the  same  angle  of  friction,  small  changes  in  cratering  efficiency 
would  be  expected.  The  results  are  consistent  with  this  expectation,  except  for 
shots  1 7- X  and  17-0.  These  two  shots,  while  consistent  with  each  other,  seem  to 
be  small  compared  to  the  others.  Figure  5.3  shows  the  measured  wet  density 
versus  moisture  content.  Note  that  the  samples  with  2.7-percent  moisture  appear 
anomalous  on  this  curve  also.  It  may  be  that  sample  preparation  variations  for 
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Table  5.1.  Summary  of  material  property  tests  on  KAFB  desert  alluvium. 
(See  Appendix  B.) 


Moisture  Content 
% 


Cohesion 

Triaxial  test  Direct  Shear 


(dyne/cni  ) 


(dyne/ctTr ) 


Angle  of  Internal 
Friction 

(degree) 


<1.0 

2.8  x 

10 

2.4  (high  pressure) 

2.0  x 

10 

3. 0-4.0 

6.5  x 

10 

4.4  (high  pressure) 

1.7  x 

10 

7.2 

3.2  x 

10 
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Fig.  5.1  Variation  of  cratering  efficiency  with  moisture  content  of  KAFB  alluvium  for  constant  it. 
(equivalent  to  70  Tons  TNT). 


CONSTANT  CONFINING  PRESSURE 
AT  CRATER  FLOOR,  pgh  =  8 


Crossplot  at  constant  value  of  pgh  (from  Fig.  4.5)  for  crater 
aspect  ratio  as  a  function  of  moisture  content  of  KAFB  alluvium. 


MOISTURE  CONTENT  (percent) 

Fig.  5.3  Bulk  density  of  KAFB  alluvium  versus  moisture  content.  (Samples  were  each  centrifuged  for  10  minutes 
at  500  G  to  achieve  uniform  stable  soil  conditions;  the  shot  number  is  shown  by  each  point.) 


the  2.7-percent-moisture  experiments  explain  the  results  of  17-X  and  17-0. 

The  variation  of  aspect  ratio  with  moisture,  as  shown  in  Fig.  5.2, 
again  shows  a  consistent  trend.  As  the  moisture  content  is  decreased,  the 
crater  becomes  flatter.  For  large  moisture  content,  the  crater  becomes  more 
nearly  hemispherical.  There  was  no  consistent  variation  in  the  measured  values 
of  angle  of  friction  for  specimens  with  moisture  in  excess  of  about  2  percent. 
The  only  consistent  variation  in  strength  for  these  samples  was  the  cohesion 
measured  in  the  direct  shear  tests,  as  shown  in  Table  5.1.  This  leads  to  the 
conclusion  that  it  is  the  decreasing  cohesion  that  causes  the  flattening  of  the 
craters  as  the  moisture  content  is  decreased.  At  the  other  extreme,  the 
7.5-percent-moisture  specimen  with  the  greatest  cohesion  has  a  bowl  shape  very 
nearly  the  same  as  the  clay. 


62 


SECTION  6 

SATURATED  SAND  RESULTS 


The  feasibility  of  performing  centrifuge  cratering  experiments  in 
saturated  sand  was  assessed.  Six  shots  were  performed  in  three  different  types 
of  sand.  The  results  are  in  qualitative  agreement,  with  the  strength  theory 
developed  above  and  because  of  their  importance  to  interpreting  the  Pacific 
Proving  Ground  (PPG)  nuclear  cratering  events,  a  discussion  is  included  here. 

Figure  6.1  shows  the  results  of  these  six  experiments  as  a  plot  of  the 
cratering  efficiency  wv  versus  the  gravity-scaled-yield  parameter  ^  A1  so 
shown  are  the  results  of  six  1-G  shots  by  Piekutowski  (unpublished  data).  The 
line  representing  the  dry  Ottawa  sand  results  is  shown  for  comparison. 

At  small  values  of  tt 2 ,  the  data  show  significantly  reduced  cratering 
efficiency  as  compared  to  the  dry  sand.  Shot  28- X  at  10  G,  shows  a  much  lower 
cratering  efficiency.  However,  for  this  crater  some  lateral  washing,  due  to 
wave  motion,  was  observed  on  the  high-speed-rnovie  coveraqe  during  the  slowdown 
of  the  centri fuqe.  This  was  consistent  with  a  reduced  crater  cross  section. 
For  that  reason,  it  is  thought  to  have  partially  filled  in  an  actual  crater 
volume  that  was  larger.  This  occurrence  is  also  suggested  by  the  shape  of  that 
crater,  which  was  much  flatter  than  all  others. 

If  shot  78-X  is  discounted  (or  considered  a  lower  limit),  the  data  are 
consistent  and  show  the  same  trend  as  that  observed  for  the  4-percent-moi sture 
alluvium.  At  low  values  of  *£,  the  cratering  efficiency  is  independent  of  ^ 
This  cohesion-dominated,  cube-root  range  corresponds  to  an  "apparent"  cohesion 
(Scott,  1963).  As  the  size  parameter  ^  increases,  the  cratering  efficiency 
crosses  over  the  dry-sand  curve,  and  the  saturated-sand  cratering  efficiency  is 
greater  than  that  for  dry  sand. 

This  general  behavior  is  entirely  consistent  with  the  strength  theory 
given  in  Section  3,  if  the  apparent  cohesion  and  the  strength  envelope  of  a 
saturated  soil  is  considered.  Seed  and  Lee  (  1967)  show  that  the  total-stress 
failure  envelope  for  a  saturated  sand  under  undrained  conditions  is  distinctly 
different  than  for  drained  conditions.  The  undrained  test,  which  corresponds  to 
dynamic  phenomena,  allows  the  pore  water  to  carry  significant  portions  of  the 
total  pressure  or  to  cavitate,  which  can  influence  the  effective  strength.  For 
dense  sands  the  net  strength  envelope  has  an  apparent  cohesion  even  though  the 
dry  sand  itself  has  none.  Furthermore,  for  sufficiently  large  confining 
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Fig.  6.1  Comparison  of  material  strength  model  for  cratering  efficiency 
with  centrifuge  experimental  results  for  saturated  sand. 


pressure,  the  angle  of  internal  friction  for  the  total  stress  envelope  is  much 
smaller  than  that  of  the  dry  sand  material  itself.  Seed  and  Lee  (1967)  present 
data  for  undrained  tests  on  saturated  Sacramento  River  sand  of  various 
porosities.  For  initially  consolidated  conditions,  the  angle  of  friction  for 
large  confining  pressure  is  on  the  order  of  7-8  degrees,  as  opposed  to  37 
degrees  for  the  dry  sand.  For  unconsolidated  initial  conditions  and  large 
confining  pressure  the  angle  is  zero. 

These  considerations,  in  conjunction  with  the  strength  theory  of 
Section  3,  indicate  that  at  small  *£  a  saturated  sand  should  behave  as  a 
cohesive  material  and  have  a  cube-root,  cohesion-dominated  regime.  As  is 
increased,  a  transition  to  a  gravity-dominated  response  is  expected.  For 
sufficiently  large  ir^,  the  effective  angle  of  friction  is  zero  or  very  small, 
and  the  cratering  efficiency  should  be  greater  than  that  for  dry  sand. 

The  values  of  the  aspect  ratio  r/h  for  the  saturated  sand  centrifuge 
shots  are  shown  versus  pgh  in  Fig.  6.2.  If  the  angle  of  friction  is  small  for 
large  values  of  then  the  stability  considerations  given  for  the  clay  results 
would  also  apply  to  ^he  saturated  sand.  That  is,  for  large  sizes,  there  would 
be  a  limiting  value  to  pgh/c  and  sufficiently  large  craters  would  have  much 
greater  aspect  ratios  (r/h).  This  trend  is,  of  course,  observed  in  the  PPG 
nuclear  craters  as  seen  in  the  comparison  plot  for  aspect  ratio  in  Fig.  6.3.  To 
facilitate  the  comparison  between  the  centrifuge-formed  craters  and  the  field 
data,  a  scaled  abscissa  is  used.  For  the  field  shots  G  equals  one  and  the 
actual  depth  is  shown.  For  the  centrifuge  shots  the  actual  crater  depth  is 
multiplied  by  test  condition  G  value.  The  intent  of  this  figure  is  to  illustrate 
how  the  crater  aspect  ratio  increases  with  crater  depth  for  the  saturated  sand. 


65 


I  ifrlHMB-MUtt  -Vi- — -  '■■'t**'" 


(  H/J )  OllVb  !D3dSV 


Fig.  6.3  Comparison  of  centrifuge  results  for  crater  aspect  ratio  with  data  from  PPG  nuclear 
and  high-explosive  craters. 


SECTION  7 
CONCLUSIONS 


The  simulation  of  large-scale  cratering  events  at  subscale  on  a 
centrifuge  has  been  successfully  extended  to  a  variety  of  important  materials. 
These  include  desert  alluvium  with  various  moisture  contents,  water-saturated 
sand,  and  an  oil-base  clay.  The  results  clearly  illustrate  the  dependence  of 
crater  volume  and  shape  on  gravity  and,  consequently,  the  need  for  elevated 
gravity  in  order  to  correctly  model  large-scale  events  at  subscale.  The 
accumulated  data,  in  conjunction  with  previous  1-G  results,  have  identified 
three  differing  mechanisms  that  dominate  in  different  size  regimes.  The  results 
have  been  unified  into  a  soi  1 -strength  mathematical  model  that  can  be  used  to 
predict  cratering  in  a  variety  of  materials,  and  with  various  high-explosive 
sources. 

It  has  proved  feasible  to  conduct  tests  in  water-saturated  sands.  At 
small-scaled  sizes,  the  cratering  efficiency  in  this  material  is  less  than  that 
in  dry  sand.  This  is  attributed  to  the  "apparent  cohesion"  of  a  dense  saturated 
granular  medium  under  undrained  conditions.  At  large  scaled-sizes,  the 
cratering  efficiency  exceeds  that  of  dry  sand.  This  is  attributed  to  the  small 
friction  angle  of  the  material  when  the  pressure  is  in  excess  of  the  critical 
confining  pressure  of  the  sand.  These  interpretations  are  consistent  with  the 
soi 1 -strength  model  that  has  been  presented. 

Tests  have  also  been  performed  in  desert  alluvium  with  various 
controlled  moisture  contents.  The  data  show  a  strength-dominated  regime  at 
small-scaled  sizes  and  become  more-or- 1  ess  asymptotic  to  the  dry-sand  behavior 
at  scaled  explosi ve- source  sizes  in  the  kiloton  and  above  range.  Small 
variations  in  the  moisture  content  had  a  significant  effect  on  both  crater  size 
and  structure,  consistent  with  the  measured  mater i al - strength  properties  and  the 
model  presented. 

Material  property  differences  also  produced  large  differences  in  crater 
shape  and  structure.  The  dry  sand  craters  had  essentially  constant  shape  with 
varying  scaled  size.  The  alluvium  craters,  at  a  fixed  moisture  content,  showed 
a  trend  toward  shallower  craters  at  decreasing  size.  For  fixed  scaled-source 
size,  shallower  craters  were  observed  with  decreasing  moisture  content. 

The  most  dramatic  change  of  shape  with  size  was  observed  for  the  clay 
and  the  saturated-sand  craters.  In  both  cases  there  was  a  threshold  above  which 
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the  craters  became  increasingly  shallow.  This  threshold  has  been  interpreted  as 
a  slope-stability  limit.  There  is  also  evidence  from  the  clay-crater  data  that 
large-scale  shallow  craters  do  not  necessarily  go  throuqh  a  transient  deeper 
bowl-shaped  phase  as  commonly  suggested. 

The  shape  of  the  saturated-sand  craters  is  consistent  with  the  PACE 
field  experiments.  Extrapol ati on  of  the  observed  reduction  of  aspect  ratio 
(radius/depth)  with  yield  is  consistent  with  the  trend  observed  for  the  PPG 
nuclear  craters. 

This  experimental  study  successfully  examined  the  appl  icabi 1 ity  of  the 
centrifuge  technique  to  a  variety  of  real-world  materials.  The  material 
properties  were  sufficiently  well  controlled  and  characteri zed  so  that  a 
meaningful  mathematical  model  that  predicts  the  effects  of  material  properties 
on  cratering  behavior  has  been  generated.  Its  validity  has  been  remonstrated 
over  a  range  of  more  than  nine  orders  of  magnitude  of  scaled  yield  and  the  model 
can  be  used  as  a  high-explosive  baseline  for  nuclear  yields  of  interest. 
Included  in  the  scaled-yield  parameter  is  a  dependence  upon  source  energy 
density.  The  validity  of  this  parameter  has  been  demonstrated  over  a  range  of 
four  to  one  in  enerqy  density  for  high  explosives.  For  impact  cratering,  a  range 
of  more  than  two  orders  of  magnitude  on  enerqy  density  is  in  agreement  with  this 
same  scaled-yield  parameter  (Schmidt  and  Holsapple,  1978a).  Furthermore,  using 
a  "working  gas"  concept  as  defined  by  Butkovich  (1967)  to  establish  a  value  for 
nuclear  energy  density  appropriate  to  the  time-scale  of  crater  excavation,  this 
source  model  promises  to  provide  a  scaling  tool  allowing  application  of  these 
centrifuge  results  to  the  prediction  of  nuclear  craters  (Schmidt  et.  aj_.  1979). 

Results  obtained  are  only  for  half-buried  charges.  They  have 
convincingly  demonstrated  that  the  centrifuqe  provides  a  viable  method  to 
determine  cratering  behavior  for  various  soils  at  large  explosive  yields  of 
interest.  In  particular  it  can  be  used  to  examine  conditions  leading  to  the 
slumping  of  transient  bowl-shape  craters.  Further,  this  is  of  particular 
interest  in  the  case  of  iturated  media  which  may  undergo  bl ast- i nduced 
liquefaction  and  thereby  reduce  the  shear  strength  to  below  that  required  for 
stabil ity. 

Recommendations  are  to  extend  the  data  base  to  include  tangent-above 
charge  configurations.  As  hei qht-of-burst  is  increased,  the  transition  between 
the  strength-dominated  cube-root  regime  and  the  gravity-dominated  regime  is 
expected  to  be  different  from  that  observed  for  half-buried  charges. 
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Preliminary  h.o.b.  data  for  scaled  yields  above  this  transition  indicate  that 
the  scaling  exponent  is  slightly  less  (1/3.8)  than  observed  for  half-buried 
explosives  (1/3.6).  This  behavior  must  be  examined  in  particular  for 
water-saturated  soils  of  varying  porosity.  Small  differences  in  the  scaling 
leads  to  large  discrepancies  in  predictions  for  large  yields.  Centrifuge 
techniques  provide  a  means  to  test  these  various  mechanisms. 
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APPENDIX  A 

CENTRIFUGE  SHOT  RECORPS 


This  section  contains  all  the  data  for  the  elevated  gravity  centrifuge 
experiments.  A  table  is  given  for  each  run  which  includes  the  test  conditions 
for  the  two  shots  performed  at  the  opposing  rotor  ends.  Following  each  table  is 
a  figure  which  contains  a  comparison  plot  for  the  two  craters  and  a  set  of 
documentary  photographs. 


A1 


SHOT  NUMBER 

OATE 

PURPOSE 

CHARGE  DESCRIPTION 
CHARGE  MASS  (qm) 

CHARGE  RADIUS  (cm) 

CHARGE  CONFIGURATION 

TEST  BED  MATERIAL 
TEST  BED  DENSITY  (qm/cc) 
MOISTURE  CONTENT  (%) 

TEST  BED  GEOMETRY 

CENTRIFUGE  SPEED  (rpm) 
GROUND  ZERO  RADIUS  (cm) 
CENTRIFUGAL  ACCELERATION 

CRATER  VOLUME  (cc) 

CRATER  RADIUS  (cm) 

MAX  CRATER  DEPTH  (cm) 
CRATER  ASPECT  RATIO  (r/h) 
CRATER  C/L  DEPTH  (cm) 

LIP  RADIUS  (cm) 

LIP  HEIGHT  (cm) 

LIP  VOLUME  (cc) 


r,V 

PIR 

PIH 


14-0 
9/19/77 
^2  Test 

CICS-1 .265(0-3) 

0.13  AqN3/l .23  PETN 
0.565 

Half-Buried  Sphere 

Permoplast  Clay  (#5 
1.53 

Homoqeneous 

571 

124.5 

(G)  454 

75.4 

4.17 

2.59 

1.61 

2.59 

5.05 

0.90 

35 

7.59E-6 

98.0 

4.46 

2.92 


1 4  -  X 
9/19/77 
*2  Test 

CICS-4(B-5) 

0.13  AqN 3/ 3 .94  PETN 
0.826 

Half-Buried  Sphere 

)  Permoplast  Clay  (#47) 
1 .53 

Homoqeneous 

471 

124.5 

309 

244 
5.91 
4.23 
1.40 
4.23 
7.20 
1 .22 

79.3 

7.26E-6 

93.4 
4.26 
3.04 
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SHOT  NUMBER  14-0  14-X 

R_ANGE  (cm)  DEPTH  (cm)  DEPTH  (cm) 


0.0 

2. 

0.6 

2. 

1.2 

2. 

1.8 

2. 

2.4 

1. 

3.0 

1. 

3.6 

0. 

4.2 

-0. 

4.8 

-0. 

5.4 

-0. 

6.0 

-0. 

6.6 

-0. 

7.2 

-0. 

7.8 

-0. 

8.4 

-  - 

9.0 

— 

9.6 

-0. 

10.2 

-  - 

10.8 

— 

11.4 

-0. 

12.0 

12.6 

13.2 

-0. 

13.8 

14.4 

— 

15.0 

-0. 

15.6 

16.2 

16.8 

-0. 

17.4 

18.0 

— 

18.6 

19.2 

—  _ 

19.8 

20.4 

_  _ 

21.0 

594 

4.231 

526 

3.986 

397 

3.942 

226 

3.826 

923 

3.575 

440 

3.196 

774 

2.948 

045 

2.366 

749 

1.590 

824 

0.800 

430 

-0.136 

046 

-0.829 

049 

-1.223 

059 

-0.859 

-0.212 

0.008 

065 

0.022 

097 

0.055 

134 

0.047 

201 

0.025 

311 

0.008 

A3 


Comparison  of  crater  profiles. 
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SHOT  NUMBER 

20-0 

20-X 

DATE 

5/18/78 

5/18/78 

PURPOSE 

1 0 -G  Check  on 

1 0 -G  Check  on 

UDRI  Data 

UDRI  Data 

CHARGE  DESCRIPTION 

C I  LAS- 1 3 ( B- 3 ) 

C ICS-5 (B-5) 

CHARGE  MASS  (qm) 

1.71  PbNfi 

0.13  AqN3/0.36  PETN 

CHARGE  RADIUS  (cm) 

0.508 

0.390 

CHARGE  CONFIGURATION 

Half-Buried  Sph  re 

Hal f-Buried  Sphere 

TEST  BED  MATERIAL 

KAFR  D. A.  #10/011 

KAFR  D. A.  #8/#9 

TEST  BED  DENSITY  (qm/cc) 

1.582 

1 .579 

MOISTURE  CONTENT  (%) 

3.9 

4.0 

TEST  BED  GEOMETRY 

Hoinoqeneous 

Homoqeneous 

CENTRIFUGE  SPETD  (rpm) 

84.5 

84.6 

GROUND  ZERO  RADIUS  (on) 

125 

1?5 

CENTRIFUGAL  ACCELERATION  (G) 

10 

10 

CRATER  VOLUME  (cc) 

43.4 

26.1 

CRATER  RADIUS  (cm) 

5.36 

4.61 

MAX  CRATER  DEPTH  (cm) 

1.37 

1 .28 

CRATER  ASPECT  RATIO  (r/h) 

3.91 

3.60 

CRATER  C/L  DEPTH  (cm) 

1.37 

1  .28 

LIP  RADIUS  (cm) 

7.80 

7.20 

LIP  HEIGHT  (cm) 

0.16 

0.14 

LIP  VOLUME  (cc) 

48.3 

42.4 

PI2 

6.08E-7 

1 .37E-7 

PIV 

40.4 

84.1 

PIR 

5.23 

6.81 

P[H 

1  .34 

1  .89 
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SHOT  NUMBER 

20-0 

20-X 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.36? 

1.283 

0.6 

1.306 

1 .251 

1.2 

1.126 

0.958 

1  .8 

0.822 

0.583 

2.4 

0.777 

0.472 

3.0 

0.653 

0.439 

3.6 

0.507 

0.29'. 

4.2 

0.?<r 

0.111 

4.8 

u.158 

-0.053 

5.4 

-0.017 

-0.105 

6.0 

-0.070 

-0.107 

6.6 

-0.1  11 

-0.119 

7.2 

-0.149 

-0,139 

7 . 8 

-0.164 

-0.133 

8.4 

-0.145 

-0.103 

9.0 

-0. i 17 

-0.100 

9.6 

-0.113 

-0.078 

’0.2 

-0.095 

-0.068 

10,  ti 

-0.081 

-0.059 

11.4 

-0.076 

-0.050 

12.0 

-0.066 

-0.046 

12.6 

-0.037 

-0.031 

13.2 

-0.036 

-i  02/ 

13.8 

-0.026 

-0.022 

14.4 

-0.021 

-0.020 

15.0 

-0.019 

-0.014 

15.6 

-0.007 

-0.014 

16  2 

-0.005 

-0.014 

7  6.8 

-0.001 

-C.002 

17.4 

0.000 

-0.002 

18. 0 

0.006 

-0.004 

18.6 

O.OOG 

-0.005 

1  °.  2 

19.8 

0.014 

-0.002 

-0.009 

20.4 

2  i  .0 

A8 

-0.018 
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Comparison  of  crater  profiles. 


Shot  20-X,  crater  formed  at  10  G  by  spherical  0.49-fjm  I’FITN 
charge’  hall  -buried  in  KAI B  desert  alluvium. 
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SHOT  NUMBER 

21-0 

2 1  -  X 

DATE 

5/26/78 

5/26/78 

PURPOSE 

Maximum  tt^ 

Maximum 

CHARGE  DESCRIPTION 

C  ICS-4(B-1 2) 

CILAS-1 3 ( B - 4 ) 

CHARGE  MASS  (qm) 

0.13  A  iN3/3 .94  PETN 

1.71  PbNg 

CHARGE  RADIUS  (cm) 

0.826 

0.508 

CHARGE  CONFIGURATION 

Half-Buried  Sphere 

Half-Buried  Sphere 

TEST  BED  MATERIAL 

KAFB  D.A.  #11 

KAFB  D.A.  #12 

TEST  BED  DENSITY  (qm/cc) 

1 .596 

1 .584 

MOISTURE  CONTENT  (%) 

4.1 

4.3 

TEST  BED  GEOMETRY 

Homoqeneous 

Homoqeneous 

CENTRIFUGE  SPEED  (rpm) 

610 

610 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G) 

520 

520 

CRATER  VOLUME  (cc) 

103 

23.7 

CRATER  RADIUS  (cm) 

5.99 

3.69 

MAX  CRATER  DEPTH  (cm) 

2.37 

1 .58 

CRATER  ASPECT  RATIO  (r/h) 

2.53 

2.34 

CRATER  C/L  DEPTH  (cm) 

2.37 

1 .58 

LIP  RADIUS  (cm) 

7.20 

4.80 

LIP  HEIGHT  (cm) 

0.32 

0.19 

LIP  VOLUME  (cc) 

47.7 

12.6 

PI2 

1 .22E-5 

3.16E-5 

PIV 

40.3 

22.1 

PIR 

4.38 

3.60 

P,H 

1.73 

1.54 
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SHOT  NUMBER 

21-0 

21- X 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.374 

1.579 

0.6 

2.351 

1.568 

1.2 

2.367 

1.282 

1.8 

2.237 

0.696 

2.4 

1 .908 

0.548 

3.0 

1.333 

0.334 

3.6 

0.945 

0.038 

4.2 

0.737 

-0.188 

4.8 

0.504 

-0.189 

5.4 

0.245 

-0.113 

6.0 

-0.002 

-0.066 

6.6 

-0.215 

-0.037 

7.2 

-0.315 

-0.020 

7.8 

-0.265 

-0.013 

8.4 

-0.183 

-0.004 

9.0 

-0.124 

-0.003 

9.6 

-0.092 

-0.002 

10.2 

-0.060 

0.001 

10.8 

-0.031 

0.004 

11.4 

-0.031 

0.011 

12.0 

-0.015 

0.012 

12.6 

-0.012 

0.014 

13.2 

-0.010 

0.01? 

13,8 

-0.006 

0.006 

14.4 

-o.oio 

0.004 

15.0 

-0.003 

0.010 

15.6 

0.00? 

0.005 

16.2 

o.nm 

0.010 

16.8 

-0.002 

0.001 

17.4 

-0.000 

0.001 

18.0 

-0.004 

0.001 

18.6 

-0.006 

0.002 

19.2 

-0.005 

0.002 

19.8 

-0.007 

0.007 

20.4 

-- 

-0.001 

21.0 

__ 
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Comparison  of  crater  profiles. 


SHOT  NUMBER 

22-X-l 

22-X-2 

DATE 

6/2/78 

6/2/78 

PURPOSE 

Stress  Gages 

Stress  Gages 

CHARGE  DESCRIPTION 

CILAS-I 3(R-5) 

CILAS-13(B-6) 

CHARGE  MASS  (qm) 

1 .70  PbN6 

1  .70  PbNfi 

CHARGE  RADIUS  (cm) 

0.508 

0.508 

CHARGE  CONFIGURATION 

Half-Buried  Sphere 

Half-Buried  Sphere 

TEST  BED  MATERIAL 

Permoplast  Clay  (#53) 

Permoplast  Clay  (#55) 

TEST  BED  DENSITY  (qm/cc) 

1 .53 

1  .53 

MOISTURE  CONTENT  (%) 

-- 

TEST  BED  GEOMETRY 

Homogeneous 

Homogeneous 

CENTRIFUGE  SPEED  (rpm) 

613 

82 

GROUND  ZERO  RADIUS  (cm) 

132.7 

132.7 

CENTRIFUGAL  ACCELERATION  (G) 

557 

10 

CRATER  VOLUME  (cc) 

40.4 

53.4 

CRATER  RADIUS  (cm) 

3.32 

3.24 

MAX  CRATER  DEPTH  (cm) 

2,37 

2.92 

CRATER  ASPECT  RATIO  (r/h) 

1.40 

1  .11 

CRATER  C/L  DEPTH  (cm) 

2.37 

2.92 

LIP  RADIUS  (cm) 

4.20 

4.20 

LIP  HEIGHT  (cm) 

0.52 

1  .09 

LIP  VOLUME  (cc) 

12.9 

51.3 

P 1  2 

3.39E-5 

6.08E-7 

PIV 

36.4 

48.1 

PIR 

3.21 

3.13 

PIII 

2.29 

2.82 
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SHOT  NUMBER 

22-X-l 

22-X-2 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.372 

2.918 

0.6 

2.283 

2.807 

1 .2 

2.042 

2.596 

1.8 

1.727 

2.266 

2.4 

1  .168 

1 .661 

3.0 

0.362 

0.583 

3.6 

-0.382 

-0.810 

4.2 

-0.523 

-1.088 

4.8 

-0.009 

-0.289 

5.4 

0.005 

-0.172 

6.0 

0.034 

-0.115 

6.6 

0.058 

-0.076 

7.2 

0.079 

-0.050 

7.8 

0.084 

-0.048 

8.4 

0.101 

-0.027 

9.0 

0.089 

-0.027 

9.6 

0.070 

-0.010 

10.2 

0.033 

0.009 

10.8 

-0.005 

0.008 

11.4 

-- 

-- 

12.0 

-- 

-- 

12.6 

-- 

-- 

13.2 

-- 

-- 

13.8 

-- 

-- 

14.4 

-- 

-- 

15.0 

-- 

-- 

15.6 

-- 

-- 

16.2 

-- 

-- 

16.8 

-- 

-- 

17.4 

-- 

-- 

18.0 

-- 

-- 

18.6 

-- 

-- 

19.2 

-- 

-- 

19.8 

-- 

20.4 

-- 

-- 

21  .0 

_ 

_  _ 

A18 


Fig.  A10.  Comparison  of  crater  profiles. 


I  uj.  Al;\  Shot  22-X-U  crater  formed  at  !»!j /  C  by  sphpric.il  1 .  /0-gm  PbN, 
charge  half- burl (?d  in  "Pennoplast"  nil-base  clay.  ° 


AX1 


Fig.  A13. 


trigger  shown  at  4  o  clock  position. 
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SHOT  NUMBER 

23-0 

23-X 

DATE 

7/14/78 

7/14/78 

PURPOSE 

Alltv 

Alluvium  Strength 

CHARGE  DESCRIPTION 

CICS-‘ji«-?» 

0 1 uS-4.0 (B-l 3) 

CHARGE  MASS  (qm) 

0.13  AgNvO.C. 

0.13  AgN3/3.95  PETN 

CHARGE  RADIUS  (cm) 

0.390 

0.826 

CHARGE  CONFIGURATION 

Hal f-Buriu  !  Sphere 

Half-Buried  Sphere 

TEST  BED  MATERIAL 

KAFB  D. A.  #1 R/#2R 

KAFB  D.A.  #3R/#4R 

TEST  BED  DENSITY  (gm/cc) 

1 .566 

1 .570 

MOISTURE  CONTENT  {%) 

3.8 

4.0 

TEST  BED  GEOMETRY 

Homogeneous 

Homogeneous 

CENTRIFUGE  SPEED  (rpm) 

358 

258 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G) 

179 

93 

CRATER  VOLUME  (cc) 

19.1 

169 

CRATER  RADIUS  (cm) 

3.96 

7.52 

MAX  CRATER  DEPTH  (cm) 

1.28 

2.50 

CRATER  ASPECT  RATIO  (r/h) 

3.09 

3.01 

CRATER  C/L  DEPTH  (cm) 

1 .22 

2.48 

LIP  RADIUS  (cm) 

4.80 

9.60 

LIP  HEIGHT  (cm) 

0.11 

0.40 

LIP  VOLUME  (cc) 

8.58 

127 

PI2 

2.45E-6 

2.18E-6 

PIV 

61.0 

65.1 

PIR 

5.83 

5.47 

PIH 

1.89 

1  ,82 
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SHOT  NUMBER 

23-0 

23-X 

RANGE  (cm) 

DEPTH  (cm) 

OEPTH  (a 

0.0 

1.221 

2.482 

0.6 

1.282 

2.495 

1 .2 

0.967 

2.503 

1.8 

0.644 

2.470 

2.4 

0.406 

2.226 

3.0 

0.223 

1.744 

3.6 

0.072 

1.415 

4.2 

-0.048 

1 .151 

4.8 

-0.107 

0.952 

5.4 

-0.101 

0.799 

6.0 

-0.063 

0.592 

6.6 

-0.039 

0.337 

7.2 

-0.018 

0.101 

7.8 

-0.006 

-0.092 

8.4 

-0.004 

-0.263 

9.0 

0.005 

-0.355 

9.6 

0.008 

-0.399 

10.2 

0.012 

-0.371 

10.8 

0.015 

-0.311 

11.4 

0.016 

-0.252 

12.0 

0.020 

-0.192 

12.6 

0.018 

13.2 

0.022 

-0.127 

13.8 

0.023 

-0.101 

14.4 

0.028 

-0.077 

15.0 

0.024 

-0.062 

15.6 

0.028 

-0.047 

16.2 

0.033 

-0.038 

16.8 

0.033 

-0.034 

17.4 

0.038 

-0.030 

18.0 

0.032 

-0.028 

18.6 

0.037 

-0.021 

19.2 

0.037 

19.8 

0.028 

20.4 

0.030 

-0.024 

Fig.  A15.  Comparison  of  crater  profiles. 


1g.  A16. 


Shot  23-0,  crater  formed  at  179  G  by  spherical  0.49-gm  PETN 
charge  half-buried  in  KAFB  desert  alluvium. 
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Fig.  A17.  Shot  23-X,  crater  formed  at  93  G  by  spherical  4,08-gm  PETN 
charge  half-burled  In  KAFB  desert  alluvium. 


SHOT  NUMBER 

24-0 

24  -  X 

DATE 

8/3/78 

8/3/78 

PURPOSE 

Effect  of  Moisture 

Effect  of  Moisture 

(Dry) 

(8%) 

CHARGE  DESCRIPTION 

C ICS-4 (R-l 5) 

C I CS-4 ( B- 1 4 ) 

CHARGE  MASS  (cjn) 

0.13  AqN3/3. 95  PETN 

0.13  AqN 3 / 3 . 96  PETN 

CHARGE  RADIUS  (cm) 

0.826 

0.826 

CHARGE  CONFIGURATION 

Half-Buried  Sphere 

Half-Buried  Sphere 

TEST  BED  MATERIAL 

KAFR  D. A.  #11/#12 

KAFB  D.A.  #7/#8 

TEST  BED  DENSITY  (gm/cc) 

1.540 

1.657 

MOISTURE  CONTENT  (%) 

<0.7 

7.5 

TEST  RED  GEOMETRY 

Homogeneous 

Homogeneous 

CENTRIFUGE  SPEED  (rpm) 

468 

468 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G) 

306 

306 

CRATER  VOLUME  (cc) 

108 

136 

CRATER  RADIUS  (cm) 

6.71 

6.10 

MAX  CRATER  DEPTH  (cm) 

1.91 

3.25 

CRATER  ASPECT  RATIO  (r/h) 

3.51 

1 .88 

CRATER  C/L  DEPTH  (cm) 

1.91 

3.25 

LIP  RADIUS  (cm) 

8.40 

7.20 

LIP  HEIGHT  (cm) 

0.39 

0.32 

LIP  VOLUME  (cc) 

90.0 

54.1 

PI2 

7.19E-6 

7.19E-6 

PIV 

40.7 

55.2 

PIR 

4.85 

4.52 

PIH 

1 .38 

2.41 
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SHOT  NUMBER 
RANGE  (cm) 
0.0 
0.6 
1.2 
1.8 

2.4 
3.0 

3.6 

4.2 

4.8 

5.4 

6.0 

6.6 

7.2 

7.8 

8.4 
9.0 

9.6 

10.2 

10.8 

11.4 

12.0 

12.6 

13.2 

13.8 

14.4 
15.0 

15.6 

16.2 

16.8 

17.4 
18.0 

18.6 
19.2 
19.8 
20. A 
21.0 


24-0 

DEPTH  (cm) 
1.912 
1.877 
1.836 
1.661 
1.439 
1.225 
1.031 
0.862 
0.672 
0.481 
0.272 
0.042 
-0.197 
-0.37? 
-0.388 
-0.304 
-0.227 
-0.165 
-0.128 
-0.090 
-0.071 
-0.051 
-0.038 
-0.049 
-0.037 
-0.038 
-0.03b 
-0.028 
-0.032 
-0.026 
-0.027 
-0.023 
-0.020 
-0.023 
-0.026 
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24- X 

DEPTH  (cm) 
3.252 
3.236 
3.147 
2.915 
2.685 
2.196 
1 .129 
0.754 
0.538 
0.300 
0.038 
-0.180 
-0.321 
-0.307 
-0.222 
-0.157 
-0.111 
-0.080 
-0.052 
-0.052 
-0.036 
-0.024 
-0.021 
-0.017 
-0.013 
-0.012 
-0.007 
-0.015 
-0.015 
-0.015 
-0.009 
-0.005 
-0.002 
-0.011 
0.001 


Fig.  A18.  Comparison  of  crater  profiles. 


Fig.  A19.  Shot  24-0,  crater  formed  at  306  G  by  spherical  4.08-gni  PETN 
charge  half-buried  in  KAFB  desert  alluvium  (< 0.7-percent 
moisture  content).  Note  flat  shallow  crater  shape  and 
compare  with  Fig.  A20. 


r 

j 

\  Fig.  A20.  Shot  24-X,  crater  formed  at  306  G  by  spherical  ^.08-gm  PETN 

:  charge  half-ouried  in  KAFB  desert  alluvium  (7.5-percent 

I  moisture  content).  Note  deep  crater  with  bench  in  wall  end 

\  compare  with  Fig.  A19. 


SHOT  NUMBER 

25-0 

2  5  -  X 

DATE 

8/14/78 

8/14/78 

PURPOSE 

Small  7T 2 

Large 

CHARGE  DESCRIPTION 

C ICS-4 ( B- 1 7 ) 

C  ICS-4 (B~  1 8 ) 

CHARGE  MASS  (gn) 

0.13  AgN3/3.97  PETN 

0.13  AgN3/3.96  PETN 

CHARGE  RADIUS  (cm) 

0.826 

0.826 

CHARGE  CONFIGURATION 

Half-Buried  Sphere 

Hal f-Buri od  Sphere 

TEST  BED  MATERIAL 

Permoplast  Clay  (#58) 

Permoplast  Clay  (#59) 

TEST  BED  DENSITY  (gm/cc) 

1.53 

1 .53 

MOISTURE  CONTENT  {% ) 

-- 

TEST  BED  GEOMETRY 

Homogeneous 

Homogeneous 

CENTRIFUGE  SPEED  (rpm) 

85 

609 

GROUND  ZERO  RADIUS  (cm) 

1 24  .6 

124.8 

CENTRIFUGAL  ACCELERATION  (G) 

10.1 

517 

CRATER  VOLUME  (cc) 

498 

150 

CRATER  RADIUS  (cm) 

6.82 

5.74 

MAX  CRATER  DEPTH  (cm) 

6.01 

2.73 

CRATER  ASPECT  RATIO  (r/h) 

1  .13 

2.10 

CRATER  C/L  DEPTH  (cm) 

6.01 

2.73 

LIP  RADIUS  (cm) 

9.00 

7.20 

LIP  HEIGHT  (cm) 

2.19 

0.65 

LIP  VOLUME  (cc) 

340 

45.5 

PI2 

2.37E-7 

1 .21 E - 5 

PIV 

187 

56.2 

PIR 

4.92 

4.14 

PIH 

4.33 

1.97 
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SHOT  NUMBER 

25-0* 

25-X* 

DATE 

10/16/78 

10/16/78 

PURPOSE 

TEST  FOR  STABILITY  OF  CRATER  SHAPE 
’  (100  SEC  AT  518  G) 

CENTRIFUGE  SPEED  (rpm) 

610 

610 

GROUND  7ER0  RADIUS  (cm) 

124.6 

124.6 

CENTRIFUGAL  ACCELERATION  (G) 

518 

518 

CRATER  VOLUME  (co) 

235 

142 

CRATER  RADIUS  (cm) 

6.77 

5.72 

MAX  CRATER  DEPTH  (cm) 

3.23 

2.53 

CRATER  ASPECT  RATIO  (r/h) 

2.10 

2.26 

CRATER  Cl L.  DEPTH  (cm) 

3.23 

2.50 

LIP  RADIUS  (cm) 

7.80 

7.20 

LIP  HEIGHT  (cm) 

0.52 

0.62 

LIP  VOLUME  (cc) 

48.5 

30.7 

PI2 

1.22E-5 

1 .22E-5 

piv 

80.1 

53.2 

4.80 

4.12 

PIH 

2.33 

1.32 
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SHOT  NUMBER  25-0  25-0* 

RANGE  (cm)  DEPTH  (cm)  DEPTH  (cm) 

0.0  6.006  3.231 


0.6 

5.954 

1  .2 

5.846 

1.8 

5.690 

2.4 

5.542 

3.0 

5.252 

3.6 

4.  °-71 

4.2 

4.382 

4.8 

3.712 

5.4 

2.850 

6.0 

1.692 

6.6 

0.400 

7.2 

-0.840 

7.8 

-1.704 

8.4 

-2.125 

9.0 

-2.185 

9.6 

-2.049 

10.2 

-0.918 

10.8 

-0.287 

11.4 

-0.069 

12.0 

-0.020 

12.6 

0.020 

13.2 

0.048 

13.8 

0.064 

14.4 

0.079 

15.0 

0.091 

15.6 

0.085 

16.2 

0.081 

16.8 

0.075 

17.4 

0.048 

18.0 

-0.007 

18.6 

-- 

19.2 

-- 

19.8 

...  - 

* After  tost  tor  stability  of 

A3  6 


3.135 
3.088 
2.997 
2.839 
2.516 
2.249 
2.033 
1 .696 
1.254 
0.708 
0.146 
-0.234 
-0.520 
-0.468 
-0.230 
0.015 
0.168 
0.120 
0.105 
0.083 
0.068 
0.059 
0.053 
0.045 
0.028 
0.013 
-0.001 
-0.028 
-0.059 
-0.110 


crater  shape:  100  sec  at  518 


SHOT  NUMBER 


25-X* 


SHOT  NUMBER 

RANGE  (cm) 

25-X 

DEPTH  (cm) 

25-X* 

DEPTH  (cm) 

0.0 

2.725 

2.496 

0.6 

2.752 

2.528 

1 .2 

2.707 

2.466 

1  .8 

2.557 

2.338 

2.4 

2.402 

2.231 

3.0 

2.155 

1 .994 

3.6 

1.832 

1.700 

4.2 

1 .383 

1 .290 

4.8 

0.899 

0.862 

5.4 

0.283 

0.264 

6.0 

-0.253 

-0.264 

6.6 

-0.577 

-0.534 

7.2 

-0.654 

-0.624 

7.8 

-0.329 

-0.191 

8.4 

-0.028 

0.072 

9.0 

0.081 

0.068 

9.6 

0.078 

0.062 

10.2 

0.068 

0.044 

10.8 

0.061 

0.032 

11.4 

0.060 

0.027 

12.0 

0.042 

0.013 

12.6 

0.049 

0.008 

13.2 

0.044 

0.010 

13.8 

0.056 

0.008 

14.4 

0.058 

0.016 

15.0 

0.066 

0.015 

15.6 

0.074 

0.020 

16.2 

0.073 

0.015 

16.8 

0.059 

0.002 

17.4 

0.031 

-0.032 

18.0 

-0.010 

-0.074 

18.6 

19.2 

-- 

19.8 

_  _ 

__ 

*After  test  for  stability  of  crater  shape:  100  sec  at  518 
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Comparison  of  crater  profiles. 
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RADIAL  POSIT  ION  (cm) 


Fig.  A22.  Comparison  of  cater  profiles. 


SHOT  NUMBER 

DATE 

PURPOSE 


26-0  26-X 

9/22,78  9/22/78 

Grain  Size  Effect  Grain  Size  Effect 

(Dense  Saturated  Sand)  (Dense  Saturated  Sand) 


CHARGE  DESCRIPTION 
CHARGE  MASS  {(jn) 
CHARGE  RADIUS  (cm) 
CHARGE  CONFIGURATION 


C  HAS- 1 3  ( B  -  7 ) 

1.70  PbNc 
0 . 508 

Hal f-Buried  Sphere 


C ILA5- 1 3 (B-8) 

1.70  PbNg 
0.508 

Half-Buried  Sphere 


i F ST  BED  MATERIAL 
TEST  BED  DENSITY  (qm/cc) 
MOISTURE  CONTENT  (%) 

TEST  RED  GEOMETRY 


Sat.  Sawinq  Sand 

2.078 

>97 

Homogeneous 


Sat.  Banding  Sand 

2.056 

>97 

Homogeneous 


CENTRIFUGE  SPEED  (rpm) 

GROUND  ZERO  RADIUS  (cm) 
CENTRIFUGAL  ACCELERATION  (G) 


265 

265 

127.1 

127.1 

100 

100 

CRATER  VOLUME  (cc) 

CRATER  RADIUS  (cm) 

MAX  CRATER  DEPTH  (cm) 
CRATER  ASPECT  RATIO  (r/h) 
CRATER  C/L  DEPTH  (cin) 

LIP  RADIUS  (cm) 

! IP  HEIGHT  (cm) 

LIP  VOLUME  (cc) 


122 

54.8 

6.88 

4.80 

2.09 

1.91 

3,29 

2.51 

2.09 

1.91 

7.80 

6.00 

0.25 

0.20 

53.8 

36.6 

PIR 

PIH 


6.08E-6 

6.08E-6 

149 

66.3 

7.36 

5.11 

2.23 

2.03 
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SHOT  NUMBER 

26-0 

26- X 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.094 

1.909 

0.6 

2.011 

1.892 

1.2 

1.781 

1.675 

1.8 

1.645 

1.421 

2.4 

1.531 

1.188 

3.0 

1.388 

0.847 

3.6 

1.232 

0.539 

4.2 

1.057 

0.284 

4.8 

0.815 

-0.007 

5.4 

0.528 

-0.190 

6.0 

0.299 

-0.202 

6.6 

0.081 

-0.162 

7.2 

-0.1 16 

-0.136 

7.8 

-0.249 

-0.112 

8.4 

-0.229 

-0.095 

9.0 

-0.191 

-0.059 

9.6 

-0  41 

-0.045 

10.2 

-0.120 

-0.031 

10.8 

-0.087 

-0.030 

11.4 

-0.066 

-0.022 

12.0 

-0.060 

-0.016 

12.6 

-0.034 

-0.015 

13.2 

-0.023 

-0.020 

13.8 

-0.021 

-0.013 

14.4 

-0.009 

-0.007 

15.0 

-0.010 

-0.006 

15.6 

-0.014 

-0.007 

16.2 

-0.008 

0.002 

16.8 

-0.006 

0.001 

17.4 

0.002 

0.004 

18.0 

0.010 

0.004 

18.6 

0.014 

0.003 

19.2 

0.024 

0.002 

19.8 

0.022 

0.004 

20.4 

0.021 

-0.001 

21.0 

0.023 

-0.008 

I 


Fig.  A29.  Shot  26-0,  crater  formed  at  100  G  by  spherical  1.70-gm  PbNc 
charge  half-buried  in  nearly-saturated  dense  Ottawa  Sawing0 
sand.  Water  has  been  drained  to  a  level  be^ow  crater  floor 
to  facilitate  crater  measurement. 
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Fig.  A30. 


sph.‘r1?a?'lC7SS™UphNh0tS9raph.°,f  crater  formed  at  100  G  by 
dens  ’S  LI;7?;^„Pb  6  ShaT  half-burled  In  nearly-satu  J. 
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Shot  26-X,  close-up  photograph  of  crater  formed  at  100  G  by 
spherical  1.70-gm  PbN^  charge  half-buried  in  nearly- saturated 
dense  Ottawa  Banding  sand.  Water  has  been  drained  to  a  level 
below  crater  floor  to  facilitate  crater  measurement. 


SHOT  NUMBER 

27-0 

27-X 

DATE 

10/9/78 

10/9/78 

PURPOSE 

Saturated  Sand 

Saturated  Sand 

Reproducibi  1  it.y 

Reproducibi 1 ity 

CHARGE  DESCRIPTION 

C I L  AS- 1 3(R-1 0) 

C I  LAS- 1 3 ( B-9 ) 

CHARGE  MASS  (qm) 

1./0  PbNfi 

1.71  PbNA 

CHARGE  RADIUS  (cm) 

0.508 

b 

0.508 

CHARGE  CONFIGURATION 

Hal f -Bur i ed  Sphere 

Hal  f-Buri ed  Sphere 

TEST  BED  MATERIAL 

Sat .  FI  i nt shot  Sand 

Sat.  Flint  shot  Sand 

TEST  BED  DENSITY  (qm/cc) 

2.113 

2.1  13 

MOISTURE  CONTENT  (%) 

>9/ 

:  97 

TEST  BED  GEOMETRY 

Homoqenoous 

"omoqeneous 

CENTRIFUGE  SPEED  (riMit) 

2  Eh 

265 

GROUND  7ERO  RADIUS  (cm) 

127.1 

127 .1 

CENTRIFUGAL  ACCELERATION  (G) 

1 00 

100 

CRATER  VOLUME  (cc, 

73.8 

86.6 

CRATER  RADIUS  (cm) 

6.03 

6.09 

MAX  CRATER  DEPTH  (cm) 

1.7/ 

1.78 

CRATER  ASPECT  RATIO  (r/h) 

3.41 

3.42 

CRATER  C/L  DEPTH  (cm) 

1.77 

1.78 

LIP  RADIUS  (cm) 

7.80 

7.80 

LIP  HEIGHT  (cm) 

0.31 

0.36 

LIP  VOLUME  (cc) 

59.2 

74.8 

PI? 

6.0HE-6 

6.08E-6 

PIV 

91.7 

108 

PIR 

6.48 

6.55 

PIH 

1  .90 

1  .91 
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SHOT  NUMBER 
RANGE  (cm) 
0.0 
0.6 
1.2 
1  .8 

2.4 
3.0 

3.6 

4.2 

4.8 

6.4 

6.0 

6.6 

7.2 

7.8 

8.4 
9.0 

9.6 
10.2 
10.8 

11.4 
12.0 
12.6 

13.2 

13.8 

14.4 
16.0 

15.6 

16.2 

16.8 

17.4 
18.0 

18.6 
19.2 
19.8 

20.4 

21 .0 


27-0 

DEPTH  (cm) 
1.766 
1.701 
1.548 
1.408 
1.249 
1 .056 
0.799 
0.544 
0.334 
0.162 
0.007 
-0.157 
-0.296 
-0.313 
-0.231 
-0.173 
-0.113 
-0.087 
-0.067 
-0.052 
-0.047 
-0.036 
-0.035 
-0.030 
-0.017 
-0.017 
-0.020 
0.003 
0.002 
0.010 
0.01? 
0.011 
0.004 
0.003 
-0.003 
-0.005 
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2  7  -  X 

DEPTH  (cm) 
1.775 
1.757 
1.644 
1 .524 
1.393 
1 .198 
0.982 
0.707 
0.437 
0.234 
0.026 
-0.161 
-0.318 
-0.363 
-0.288 
-0.225 
-0.180 
-0.123 
-0.108 
-0.095 
-0.063 
-0.038 
-0.030 
-0.027 
-0.014 
-0.015 
-0.012 
-0.005 
0.014 
0.015 
0.009 
-0.001 
0.006 
0.001 
0.000 
-0.010 


Fig.  A34.  Comparison  of  crater  profiles. 
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Fig.  A36.  Shot  27-0,  crater  formed  at  TOO  G  by  spherical  1.70-gm  PbN6 
charge  half-buried  in  nearly-saturated  dense  Ottawa 
Flintshot  sand.  At  shot  time,  water  table  was  tangent 
to  surface  at  the  charge  location. 


Fig.  A38. 


Pre-shot  photograph  (27-X)  showing  charge  placement  and 
initial  water  level. 
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Fig.  A39.  Shot  27-X,  crater  formed  at  100  G  by  spherical  1.70-gm  PbNfi 
charge  half-buried  in  nearly-saturated  dense  Ottawa  ' 

Flintshot  sand.  At  shot  time,  water  table  was  tangent 
to  surface  at  the  charge  location. 
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Fig,  A40,  This  photograph  shows  crater  ( 2 7 - X )  after  water  was  drained 
to  a  level  below  crater  floor  to  facilitate  measurement. 

See  Fig.  A39. 
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SHOT  NUMBER 

28-0 

Z8-X* 

DATE 

11/2/78 

11/2/78 

PURPOSE 

Hi gh-G 

Low-G 

CHARGE  DESCRIPTION 

CILAS(B-H) 

C I LAS ( B- 1  3) 

CHARGE  MASS  (gm) 

1.70  PbNfi 

1.70  PbNg 

CHARGE  RADIUS  (cm) 

0.508 

0.508 

CHARGF  CONFIGURATION 

Half-Buried  Sphere 

Half-Buried  Spher& 

TEST  RED  MATERIAL 

Sat.  Flintshot  Sand 

Sat.  Flintshot  Sand 

TEST  BED  DENSITY  (qtr/cc) 

2.107 

2.109 

MOISTURE  CONTENT  (%) 

>97 

>97 

TEST  BED  GEOMETRY 

Homogeneous 

Homogeneous 

CENTRIFUGE  SPEED  (rpm) 

605 

84 

GROUND  ZERO  RADIUS  (cm) 

127.1 

127.1 

CFNTRIFUGAL  ACCELERATION  (G) 

500 

10 

CRATER  VOLUME  (cc) 

34.1 

67.5 

CRATER  RADIUS  (cm) 

4.68 

5.81 

MAX  CRATER  DEPTH  (cm) 

1  .13 

1  .45 

CRATER  ASPECT  RATIO  (r/h) 

4.14 

4.01 

CRATER  C/L  DEPTH  (cm) 

1  .13 

1  .30 

LIP  RADIUS  (cm) 

6.00 

7.80 

LIP  HEIGHT  (cm) 

0.27 

0.12 

LIP  VOLUME  (cc) 

36.2 

15.5 

pi2 

3.16E-5 

6.08E-7 

f'v 

42.3 

83.7 

-R 

5.03 

6.24 

PV 

1.21 

1.56 

*'l  diameter  only 
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SHOT  NUMBER 

28-0 

28-X* 

RANGE  (cm) 

DEPTh  (cm) 

DEPTH  (cm) 

0.0 

1.129 

1.299 

0.6 

1.018 

1 .455 

1.2 

0.862 

1.445 

1 .8 

0.779 

1.339 

2.4 

0.698 

1.176 

3.0 

0.595 

0.895 

3.6 

0.453 

0.722 

4.2 

0.197 

0.510 

4.8 

-0.050 

0.264 

5.4 

-0.198 

0.124 

6.0 

-0.270 

-0.014 

6.6 

-0.206 

-0.052 

7.2 

-0.138 

-0.072 

7.8 

-0.118 

-0.124 

8.4 

-0.0/8 

-0.061 

9.0 

-0.070 

-0.022 

9.6 

-0.055 

-0.031 

10.2 

-0.039 

-0.017 

10.8 

-0.035 

-0.046 

11.4 

-0.018 

-0.022 

12. C 

-0.009 

-0.021 

12.6 

-0.005 

0.022 

13.2 

-0.006 

0.016 

13.8 

-0.008 

0.013 

14.4 

-0.001 

0.038 

15.0 

0.002 

0.035 

15.6 

0.002 

0.030 

16.2 

0.002 

0.025 

16.8 

0.002 

0.033 

17.4 

0.004 

0.030 

18.0 

0.002 

0.030 

10.6 

0.003 

0.017 

19.2 

0.003 

0.024 

19.8 

0.001 

0.013 

20.4 

-0.002 

0.008 

21.  U 

-0.020 

-0.008 
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Fig.  A41 .  Comparison  of  crater  profiles. 


rig.  A42. 


Shot  28-0,  crater  formed  af  con  r  k 

charge  half-burled  In  ne»r'VSatura?JP i"' Hcal  PbN, 

Hintshot  sand.  y  saturated  dense  Ottawa  6 


A65 


I  Icj.  A44. 


Shot  2H-X,  crater  formed  at  10  0  by  spherical  1.70- qm  PbN^ 
charge  half-hurled  In  nearly-saturated  dense1  Ottawa 
I'Untshot  sand.  At  shot  time,  water  table  was  tangent 
to  surface  at  the  charge  location. 
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Shot  2H-X,  close-up  photograph  showing  effects  of  water  wave 
washing  of  crater  walls.  Z-axis  shown  on  photo  Is  line  of 
tangency  of  water  table  In  flight.  See  I  1 g .  A44, 
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SHOT  NUMBER 

29-0 

29  -  X 

DATE 

1 1/16/78 

11/16/78 

PURPOSE 

Ground  Motion 

Ground  Motion 

CHARGE  DESCRIPTION 

C  ICS-4 (B-20) 

CICS-5(B-6) 

CHARGE  MASS  (qm) 

0.13  AqNj/3.96  PETN 

0.13  AqN 3 /0 .36  PETN 

CHARGE  RADIUS  (cm) 

0.826 

0.390 

CHARGE  CONFIGURATION 

Half-Buried  Sphere 

Half-Buried  Sphere 

TEST  RED  MATERIAL 

Permoplast  Clay  (#61) 

Permoplast  Clay  (#60) 

TEST  RED  DENSITY  (qm/cc) 

1 .63 

1 .53 

MOISTURE  CONTENT  {%) 

-- 

TEST  RED  GEOMETRY 

Homoqeneous 

Homoqeneous 

CENTRIFUGE  SPEED  (rpm) 

460 

602 

GROUND  ZERO  RADIUS  (cm) 

124.5 

124.5 

CENTRIFUGAL  ACCELERATION  (G) 

294 

504 

CRATER  VOLUME  (cc) 

211 

17.6 

CRATER  RADIUS  (cm) 

6.70 

2.56 

MAX  CRATER  DEPTH  (cm) 

3.7/ 

1  .84 

CRATER  ASPECT  RATIO  (r/h) 

1.51 

1 .39 

CRATER  C/L  DEPTH  (cm) 

3.77 

1  .84 

LIP  RADIUS  (cm) 

7.20 

3.00 

LIP  III  I GILT  (cm) 

0.90 

0.40 

LIP  VOLUME  (cc) 

7  7.9 

9.93 

PI, 

6.90E-6 

6. 9 HE -6 

P!V 

79.1 

55.0 

pi 

4.11 

3.74 

(,I|| 

2.72 

2.69 
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SHOT  NUMBER 
RANGE  (cm) 
0.0 
0.6 
1  .2 
1.8 

2.4 
3.0 

3.6 

4.2 

4.8 

5.4 

6.0 

6.6 

7.2 

7.8 

8.4 
9.0 
9.6 

10.2 

10.8 

11.4 
12.0 
12.6 

13.2 

13.8 

14.4 
15.0 

15.6 

16.2 

16.8 

17.4 
18.0 

18.6 
19.2 
19.8 

20.4 

21.0 


29-0 

DEPTH  (cm) 
3.767 
3.772 
3.696 
3.548 
3.320 
2.956 
2.527 
2.029 
1.300 
0.392 
-0.385 
-0.752 
-0.903 
-0.589 
-0.285 
0.014 
0.050 
0.065 
0.070 
0.083 
0.079 
0.078 
0.067 
0.046 
0.041 
0.020 
-0.001 
-0.022 
-0.054 
-0.095 
-0.150 


29  -  X 

DEPTH  (cm) 
1 .839 
1.638 
1 .358 
0.884 
0.135 
-0.400 
-0.369 
-0.048 
0.012 
0.045 
0.071 
0.079 
0.096 
0.103 
0.119 
0.116 
0.119 
0.122 
0.119 
0.1 19 
0.110 
0.104 
0.080 
0.063 
0.046 
0.014 
-0.008 
-0.039 
-0.096 
-0.137 
-0.202 
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Fig.  A46.  Comparison  of  crater  profiles. 


Fig.  A47.  Shot  29-0,  crater  formed  at  294  G  by  spherical  4.08-gm  PETN 
charge  half-buried  in  "Permoplast"  oil-base  clay.  Gridwork, 
holes  and  sand  columns  were  used  to  record  residual  ground 
displacement  in  near  field.  Compare  with  lf2-scale  event 
shown  in  Fig.  A48. 


Fig.  A48.  Shot  29-X,  crater  formed  at  504  G  by  spherical  0.49-gm  PETN 
charge  half-buried  in  "Permopl ast"  oil-base  clay.  Gridwork, 
holes  and  sand  columns  were  used  to  record  residual  ground 
displacement  in  near  field.  Compare  with  twice-scale  event 
shown  in  Fig.  A47. 
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MATERIAL  PROPERTIES  OF  TEST  SOILS 


REPORT  ON 
LABORATORY  TESTS 
CENTRIFUGE  CRATERING  STUDY  I] 
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BOEING  AEROSPACE  COMPANY 


NOVEMBER,  1978 


SHANNON  A  WILSON,  INC. 
Geotechnical  Consultants 
1105  North  38th  Street 
Seattle,  Washington  98103 
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REPORT  ON 
LABORATORY  TESTS 
PERFORMED  ON  SOIL  MATERIALS 
FOR 

CENTRIFUGE  CRATERING  STUDY  II 


I.  PURPOSE  AND  SCOPE 


This  report  presents  the  results  of  a  series  of  unconsol idated- 
und rained  (UU)  triaxial  compression  tests  and  direct  shear  tests  performed  on 
five  soil  samples  provided  by  Boeing  Aerospace  Company.  Three  UU  tests  and  one 
direct  shear  test  were  performed  on  each  material  type  to  determine  the  shear 
strength  parameters  of  the  materials.  The  samples  were  designated  as  follows: 
Ottawa  Sand  No.  20,  Ottawa  Sand  No.  SO,  Ottawa  Sand  No.  70,  KAFBDA  and  Modeling 
Clay. 


II.  test  PROCEDURES 

A.  Unconsol idatod-Undrai ned  T^i axial  Tests 

Three  ()JU)  test  specimens  were  prepared  for  each  sample.  Test  speci¬ 
mens  of  the  sand  were  prepared  by  sifting  the  dry  sand  through  a  No.  16  mesh 
sieve  from  a  height  of  32  inches  into  a  split:  mold  to  obtain  the  density 
specified  by  your  office.  For  the  alluvium  type  soil,  KAFBDA,  the  density  of 
the  test  specimens  were  obtained  by  vibrating  a  known  weight  of  material  at  it.s 
natural  water  content  into  a  split  mold  of  known  volume.  Test  specimen  of  the 
modeling  clay  w ere  prepared  by  placing  a  known  weight,  of  clay  in  ten  equal  layers 
into  a  split  mold  and  hand  tamping  each  layer  with  a  rod  to  obtain  the  desired 
density.  The  mold  used  was  lined  with  a  thin  rubber  membrane  that  was  expanded 
against  the  split  mold  by  vacuum. 

After  each  UU  test  specimen  was  prepared,  it  was  mounted  -in  a  triaxial 
test  chamber  and  subjected  to  a  specified  triaxial  confining  pressure.  It  was 
then  sheared  under  strai n-control led  conditions,  while  maintaining  a  constant 
con+ining  pressure  and  without  allowing  drainage.  For  each  material  type, 
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triaxial  confining  pressures  of  10,  20  and  30  psi  were  used  on  the  first,  second 
and  third  test  specimens,  respectively,  os  specified.  Water  was  used  as  the 
confining  medium  Tor  all  tests. 

Simultaneous  readings  of  load  and  deformation  were  obtained  at  "egwlar 
time  intervals  throughout  the  shearing  period.  Deviator  stress  and  axial  strain 
values  were  then  computed  from  these  readings. 

B.  Direct  Shear  Tests 

Direct  shear  tes^s  were  performed  in  general  accordance  with  ASTM 
D-3080-/2.  One  direct  shear  test  was  performed  on  each  material  type.  Test 
specimens  of  the  sand  were  prepared  in  a  similar  manner  as  those  for  the  IIU 
specimens  except  they  were  prepared  in  a  4  x  4-inch  shear  box.  The  alluvium  and 
modeling  clay  samples  were  packed  in  the  box  to  the  designated  density. 

After  each  specimen  was  compacted  in  the  shear  box,  a  very  normal  load 
was  applied.  Each  specimen  was  sheared  horizontaly  in  a  constant  strain  shear 
device  under  drained  conditions. 
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A. 93 
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0.  0 

0  .  n 

0.0 

0 . 0 
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0.0 
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Fig.  B 2 

Ottawa  sand  "FLINTSH0T" 
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National  Sec  (He  Mitigation  4  Rsch 

ATTN:  Assistant  Associated  Dir 
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Lawrence  Livermore  National  Lab 
ATTN:  W.  Crowley 

ATTN:  Technical  Info  Dept,  Library 
ATTN:  L-10,  H.  Kruger 

Los  Alamos  National  Laboratory 
ATTN:  M,  Henderson 
ATTN:  MS  670,  J.  Hopkins 
ATTN:  MS362,  Llbrarlar 
ATTN:  R.  Brldwell 
ATTN:  G,  Spillman 
ATTN:  MS  364  (Class  Reports  Lib) 

Oak  Ridge  National  Laboratory 

ATTN:  Central  Rsch  Library 
ATTN:  Civil  Def  Res  Proj 

Sandla  National  Laboratories 
Livermore  Laboratory 

ATTN:  Library  &  Security  Classification  Dlv 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

California  Research  &  Technology,  Inc 
ATTN:  Library 
ATTN:  K.  Kreyenhagen 
ATTN:  S,  Schuster 
ATTN:  M.  Rosenblatt 

University  of  California 
ATTN:  J.  Cheney 

California  Research  S  Technology,  Inc 
ATTN:  D.  Orphal 

Calspan  Corp 

ATTN:  Library 

University  of  Denver 

ATTN:  J.  Wlsotskl 

EG&G  Wash.  Analytical  Svcs  Ctr,  Inc 
ATTN:  Library 


Sandla  National  Laboratories 


ATTN: 

A. 

Chabla 

ATTN: 

3141 

ATTN: 

R. 

Schmidt  5732 

ATTN: 

T. 

Bergs tresser  5533 

ATTN: 

H. 

Sutherland  4732 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 
Aerospace  Corp 

ATTN:  Technical  Information  Services 

Agbablan  Associates 

ATTN:  M.  Agbablan 

Applied  Research  Associates,  Inc 
ATTN:  J.  Bratton 
ATTN:  N.  Higgins 

Applied  Theory,  Inc 
2  cy  ATTN:  J.  Trullo 

AVCO  Research  &  Systems  Group 
ATTN:  Library  A830 

University  of  Arizona 

ATTN:  W.  McKinnon 

University  of  Arizona  State 
ATTN:  R.  Greeley 
ATTN:  J.  Fink 


Electromech  Sys  of  New  Mexico,  Inc 
ATTN:  L.  Piper 

Electromech  Sys  of  New  Mexico,  Inc 
ATTN:  R.  Shunk 

Eric  H.  Wang 

Civil  Engineering  Rsch  Fac 
University  of  New  Mexico 
ATTN:  N.  Baum 


Gard,  Inc 

ATTN:  G.  Neldhardt 

Geocenters,  Inc 

ATTN:  E.  Marram 
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ATTN:  B.  Hartenbaum 

Horizons  Technology,  Inc 
ATTN:  R.  Kruger 

I  IT  Research  Institute 

ATTN:  Documents  Library 

Institute  for  Defense  Analyses 
ATTN:  Classified  Library 

JAYCOR 

ATTN:  H.  Llnnerud 


BDM  Corp 

ATTN:  Corporate  Library 
ATTN:  T.  Neighbors 

Boeing  Co 

ATTN:  Aerospace  Library 


Boeing  Co 

4  cy  ATTN:  K,  Holsapple 
20  cy  ATTN:  R.  Schmidt 


California  Institute  of  Technology 
ATTN:  T.  Ahrens 


Kaman  AvIDyne 

ATTN:  Library 

Kaman  Sciences  Corp 
ATTN:  Library 

Kaman  Tempo 

ATTN:  DASIAC 


Lockheed  Missiles  &  Space  Co,  Inc 

ATTN:  Technical  Information  Center 

ATTN:  T,  Geers 

ATTN:  Technical  Library 
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Lockheed  Missiles  4  Space  Co,  Inc 
ATTN:  TIC-Llbrary 

Martin  Marietta  Corp 
ATTN:  G.  Freyer 

Lunar  and  Planetary  Institute 
ATTN:  P.  Schultz 

McDonnell  Douglas  Corp 
ATTN:  R.  Halprln 

Merritt  CASES,  Inc 
ATTN:  Library 
ATTN:  J.  Merritt 

University  of  New  Mexico 

ATTN:  CERF,  G.  Leigh 
ATTN:  CERF,  N.  Baum 

Pacific-Sierra  Research  Corp 
ATTN:  H.  Brode 

Pacifica  Technology 
ATTN:  G.  Kent 
ATTN:  Tech  Library 

Patel  Enterprises,  Inc 
ATTN:  M.  Patel 

Physics  Applications,  Inc 
ATTN:  F,  Ford 

Physics  International  Co 
ATTN:  J.  Thomsen 
ATTN:  F,  Sauer 
ATTN:  Technical  Library 
ATTN:  L.  Behrmann 
ATTN:  E.  Moore 

R  &  D  Associates 

ATTN:  J .  Lewis 

ATTN:  Technical  Information  Center 

ATTN:  R.  Port 

ATTN:  J.  Carpenter 

ATTN:  W.  Wright 

ATTN:  P.  Haas 

Science  Appl  Icatlons,  Inc 
ATTN:  J .  Dlshon 

Science  Applications,  Inc 

ATTN:  Technical  Library 
ATTN:  M.  Wilson 

Science  Applications,  Inc 
ATTN:  D,  Bernstein 
ATTN:  D,  Maxwell 
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Science  Applications,  Inc 
ATTN:  W.  Layson 
ATTN:  J.  Cockayne 
ATTN:  M.  Knasel 

Southwest  Research  Institute 
ATTN:  A,  Wenzel 
ATTN:  W.  Baker 

Science  4  Engineering  Associates,  Inc 
ATTN:  J.  Stockton 

SRI  International 

ATTN:  B.  Gasten 
ATTN:  G,  Abrahamson 

Systems,  Science  4  Software,  Inc 
ATTN:  T.  Rlney 
ATTN:  D.  Grlne 
ATTN:  T.  Cherry 
ATTN:  K.  Pyatt 
ATTN:  R.  Lafrenz 
ATTN:  Library 

Terra  Tek,  Inc 

ATTN:  Library 
ATTN:  S.  Green 

Tetra  Tech,  Inc 

ATTN:  L.  Hwang 

TRW  Defense  4  Space  Sys  Group 

ATTN:  Technical  Information  Center 
ATTN:  I.  Alber 
ATTN:  R.  Plebuch 
ATTN:  D.  Baer 
ATTN:  J.  O'Keefe 
?  cy  ATTN:  N.  Llpner 

TRW  Defense  4  Space  Sys  Group 
ATTN:  E.  Wong 
ATTN:  P.  Dal 

Universal  Analytics,  Inc 
ATTN:  l .  Field 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  Wright 
ATTN:  M.  Baron 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  Isenberg 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  A.  Misovec 
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